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Abstract 

4-Bromomethyl coumanns (II) have been reacted with different hydroxy coumanns to obtain 
biscoumarinyl ethers as analogues of Lasiocephalin. All the compounds have been 
characterized by spectral data and anti-microbial activity of some of them has been reported. 

(Keywords . spectral data/bis coumarinyl ethers/hydroxy coumarins/lasiocephalin). 

Introduction 

Compounds containing two coumarin moieties have been found to be extremely 
useful as anticoagulants*, antimicrobials^ and triplet sensitizers^. Further Bhattacharya 
and Das'* have isolated a biscoumarinyl ether Lasiocephalin^ (I) from the roots of 
“Lasiocepkon eriocephalns” the structure of which was deduced by Vishwanathan 
and Balakrishnan^. An Important feature of this compound is the presence of an aryl 
ether linkage between C^ and Cs position of the two coumarin moieties. In the present 
paper 6-hydroxy-4-methyl coumarin’, 4-hydroxy coumarin* and 7, 8-dihydroxy-4- 
methyl coumarin®, have been used as a nucleophiles and the ether linkage has been 
established by substituting the bromine situated at the allylic position with respect to 
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the biogenetically important C3-C4 double bond’® in various 4 -bromomethyl 
coumarins”’ 



(I) 



R = -6-CH3, -7-CH3, - 5 , 6 -Benzo, - 7 , 8 -Benzo, -7-OCH3, -6-OCH3, 6-Cl 

Scheme - 1 
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Materials and Method 

The methodology adopted for the synthesis of the title compounds is outlined in 
Scheme-1. 4-Bromomethyl coumarins'^ (II) which are the required synthons have 
been synthesised by the Pechmann cyclisation of various phenols with 4- 
bromoethylacetoacetate''* The reaction of 6-hydroxy-4-methyl coumarin with (O), 
resulted in the formation of ethers (III). However, in the reaction of 4-hydroxy 
coumarin with (II) an equimolar quantity of AgNOa was used to bring about the ether 
formation. The reaction probably proceeds via the formation of the silver salt of 4- 
hydroxy coumarin, which is evidenced by the formation of silver bromide which 
precipitates in the reaction mixture. 

Out of the two phenolic hydroxyls in 4-methyl daphnetin the C7-OH is more 
nucleophilic than Cs-OH and is also less sterically hindered. The reaction of (11) with 
this substrate has shown that the ether linkage is established with the C7-OH and the 
resulting ethers have the structure (V) with a free Cg-OH group. This is further 
supported by the formation of 7-benzyloxy-8-hydroxy-4-methyl coumarin in the 
reaction of benzyl chloride with 7,8-dihydroxy-4-methyl coumarin. Spectral and 
physical data (m.p. 180-81°) of this compound agrees with the literature report*^. 

Melting points reported are uncorrected. IR spectra (v^ax in cm’’) were recorded in 
KBr on Nicolet. Impact-410. FT-IR spectrometer. PMR spectra were recorded on 
Brucker-Varian-300 MHz FT-NMR using TMS as internal standard. The elemental 
analyses was recorded on Thermo Quest EA 1110 CHN rapid analyser. 

4-Bromoethylacetoacetate was prepared according to the method of Burger and 
Ullyot''* and substituted 4-bromomethyl coumarins (11) were synthesised by the 
Pechmann cyclisation of various phenols with 4-bromoethylacetoacetate according to 
the literature methods.'^ Amongst the hydroxy coumarins, 4-hydroxy coumarin was a 
commercial sample and 6-hydroxy-4-methyl and 7,8-dihydroxy-4-methyl coumarins 
were prepared by the reaction of hydroquinone and pyrogallol with ethyl acetoacetate 
using 70% sulfuric acid. 

I. Synthesis of 6-(0-4'-coumarino methyl)-4-methyl-coumarins (III) 

General method : 

A mixture of 6-hydroxy-4-methyl coumarin (0.01 mol), anhydrous potassium 
carbonate (0.01 mol.) was stirred for 30 min. in dry acetone (30 ml). To this, 4- 
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bromomethyl coumarin (11) (0.01 mol) was added and the stirring was continued for 
24 h. The separated solid, was filtered, and washed with 1:1 HCl (30 ml) and with 
water and recrystallised from dioxane. The compounds prepared are listed in Table 1. 


Table 1- Physical and analytical data for biscoumarinyl ethers 


No 

Compound 

M P°c 

Yield % 

Mol. formula 

Analysis(%) 

Req (found) 

R 

No 

C 

H 

! 

6 -CH 3 

Ilia 

272-73 

72 

QiHjeOs 

72 40 

4 63 







(72 20) 

(4 50) 

2 . 

7 -CH 3 

Illb 

282-84 

64 

C21H16OS 

72 40 

4 63 







(72 16) 

(4 38) 

3 

5,6- 

IIIc 

266 

78 

C24H16O5 

75 00 

4 16 


Benzo 





(74 80) 

(4 08) 

4. 

7 , 8 - 

Hid 

286-88 

80 

C24H16O5 

75 00 

4 ^6 


Benzo 





(74 78) 

(4 02) 

5 

6-OCH3 

Ille 

270 

76 

QiHieOe 

69 23 

4 39 







(69 02) 

(4.08) 

6 

7 -OCH 3 

Illf 

292-94 

70 

QiHifiOe 

69 23 

4 39 







(68 98) 

(4 18) 

7. 

6 -Cl 

Illg 

275 

72 

C 20 H 13 O 5 CI 

65 13 

3 53 







(64 96) 

(3 41) 

8 . 

6 -CH 3 

IVa 

110 

70 

C20H14O5 

71 85 

4 19 







(71 50) 

(4.10) 

9. 

5,6- 

Ivb 

160 

75 

C 23 H 14 O 5 

74 59 

3.78 


Benzo 





(74 21) 

(3.56) 

10 . 

7,8- 

Ivc 

154-55 

72 

C23H14O5 

74.59 

3 78 


Benzo 





(73.80) 

(3.60) 

11 . 

6-OCH3 

Ivd 

138 

75 

C2oH}406 

71.85 

4.19 







(71 50) 

(4 10) 


Table Contd.. 
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Table 1 Contd. 


12 

6 -Cl 

IVe 

104 

70 

C19H11O5CI 

64 32 

3 10 







(64 10) 

(3 02) 

13 

6-CH3 

Va 

250® 

80 

C21H16O6 

69 23 

4 40 







(69 10) 

(4 30) 

14. 

7 -CH 3 

Vb 

258 

78 

QiHieOe 

69 23 

4 40 







(69 05) 

(4 20) 

15. 

5.6- 

Vc 

265 

80 

C24H16O6 

72 00 

4.00 


Benzo 





(71 80) 

(3.90) 

16 

7,8- 

Vd 

278 

75 

C24H16O6 

72.00 

4 00 


Benzo 





(71 75) 

(3 88 ) 

17 

6 -OCH 3 

Ve 

232 

80 

C21H16O7 

66 32 

4.21 







(66 10 ) 

(410) 

18. 

7 -OCH 3 

Vf 

246 

70 

C21H16O7 

66 32 

421 







(66 20 ) 

(4 15) 

19. 

6 -Cl 

Vg 

262 

75 

C20H13O6CI 

62 42 

3 38 







(62 25) 

(3 18) 


Etiiers III & V were recrystallised from dioxane and IV were recrystallised from ethanol/aq ethanol 


II. Synthesis of 4-(0-4'-coumarino methyl)-coumarins (IV) 

General method : 

A mixture of 4-hydroxy coumarin (0.01 mol), and silver nitrate (0.01 mol) in 
rectified spirit (30 ml) was heated on water bath for half an hour using a reflux 
condenser. To this was added, 4-bromomethyl coumarin (H) (0.01 mol) and heated on 
a water bath for 3 h. The pale yellow precipitate of silver bromide separates out. From 
the hot solution AgBr is separated by filtration. The filtrate was added to ice cold 
water. The separated solid was washed with water, alcohol, and was recrystallised 
from ethanol or aqueous ethanol. The compounds prepared are listed in Table 1 . 
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III. Synthesis of 7-(0-4'-coumarmo methyl)-8-hydroxy-4-methyl-coumarms (V) 
General method: 

A mixture of 7,8-dihydroxy-4-methyl coumarin (0.01 mol) and anhydrous 
potassium carbonate (0.01 mol) in dry acetone (30 ml) was stirred for 1 h. To this was 
added 4-bromomethyl coumarin (H) (0.01 mol) and the stirring was contmiied for 
24 h- The separated solid was filtered and washed with 1:1 HCl and water. It was 
then dried and recrystallised from dioxane. The compounds prepared are listed in 
Table 1. 


Results and Discussion 

The biscoumarinyl ether (Ilia, R=6-CH3) obtained by the reaction of 4- 
bromomethyl coumarins and 6-hydroxy-4-methyl coumarin exhibited C=C-H 
vibrations at 2919 cm'* and the characteristic band at 1702 cm'* due to the C=0, 
stretching vibrations of coumarin. The C^C and C-O-C vibrations were observed at 
1579 cm'* and 1171 cm'* respectively (Table 2). In the PMR spectra the, methyl and 
methylene protons appeared as singlets at 2.46 5 ppm and 5.29 6 ppm respectively. 
The aromatic multiplet was observed in the region of 6.95-7.5 5 ppm. The two C3 
protons appeared as singlets, at 6.69 and 6.36 5 ppm. The C3-H2 proton which is in 
close spatial proximity with the ether oxygen appears down field when compared with 
the other C3-H1 (Table 2) which is in agreement with an earlier observation on 4- 
aryloxymethyl coumarins*^. 

Ether (IVa, R = 6-CH3) exhibited C=C-H vibrations at 2919 -cm"* and the 
characteristic band at 1715 cm'* due to C =0 stretching vibration of coumarin. The 
C=C and C-O-C vibrations were observed at 1574 cm'* and 1196 cm'* respectively. 
(Table 2). In the PMR spectra, the methyl and methylene protons appeared as singlets 
at 2.46 8 ppm and 4.68 5 ppm respectively. The aromatic multiplet was observed in 
the region of 7.28-7.43 5 ppm. The two C3 protons appeared as singlets at 6.57 5 ppm 
and 5.62 8 ppm. The C3-H2 proton which is in close spatial proximity with the ether 
oxygen appears down field when compared with the other C3-Hi (Table 2) which is in 
agreement with our earlier observations, on 4-aiyloxymethyl coumarins*^. 

Ether (Va, R = 6-CH3) exhibited a broad band at 3413 cm'* due to -OH stretching, 
C=C-H vibrations at 2919 cm'* and characteristic band at 1721 cm'* due to C =0 
stretching vibrations of coumarin. The C=C and C-O-C observed at 1572 cm'* and 
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1110 cm’’, and respectively (Table 2). In the PMR spectra, the methyl and methylene 
protons appeared as singlets at 2.38 and 2.45 5 ppm and 5.20 5 ppm respectively. The 
aromatic multiplet was observed in the region of 6.90 - 7.50 5 ppm, where as the 
phenolic -OH appeared at 5.40 5 ppm (Table 2). 

Mass Spectra : The mass spectra by electron ionisation (El) have been recorded 
for ethers IIIc, IVc and Va. The molecular ion peaks of moderate intensities have been 
observed for Va and IVc at m/z values 364 (17%) and 384 (81%) respectively. 
However, in the case of HIc the NT ion at m/z (370) is not observed. The main 
fragmentation pattern appears to be the a cleavage of ethers resulting in the formation 
of even electron ions which further follow the established fragmentation pathway of 
coumarins and benzofurans’’. A plausible fragmentation pattern is out lined for the 
benzo ether (HIc) in Scheme 2. 

In the case of the hydroxy ether Va the low intensity molecular ion at m/z 364 
17% is observed. Though the expected a-cleavage is observed in this compound, 
another predominant path way is the transfer of the phenolic hydrogen to the 
methylene group via a cyclic six membered ring leading to the formation of 4,6- 
dimethyl coumarin fragment at m/z 174 in high intensity. This odd electron ion can 
easily expel a molecule of carbon monoxide and a methyl radical leading to the 
formation of benzofuran fragments. 

Antimicrobial studies ■ Biscoumarinyl ethers III, IV and V synthesized during the 
present investigation have been subjected to a preliminary antimicrobial screening. 
The compounds were chosen at random by selecting two compounds from each HI, IV 
and V and have been screened against one bacterial and one fungal species at a dose 
of 100 pg in 10 ml of DMF as a solvent. The methodology adopted for the evaluation 
is described as follows. 

Antibacterial activity : Antibacterial activity of test compounds were evaluated 
against Pseudomonas chlororaphis bacteria by cup plate method'*. The bacteria were 
subcultured in a medium containing peptone (0.5%), yeast extract (0.15%), beef 
extract (0.15%), sodium chloride (0.35%), potassium dihydrogen phosphate (0.13%). 
Nutrient agar which served as the basal medium was prepared by dissolving 
bacteriological peptone (0.61%), yeast extract (0.3%), beef extract (0.13%) and agar 
(2.1%) in distilled wat^. The solution was sterilised for 20 min. at 15 lbs pressure in 
an autoclave. 



Table 2- IR and PMR spectral data of biscoumarinyl ethers 
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The basal medium 25-30 ml with glucose solution [(0.1 ml, 10%) to hasten the 
bacterial growth] with bacterial culture was poured in sterile petridishes. After the 
solidification of the medium holes of 9 mm diameter were bored to form cups with the 
help of a sterile cork borer. To this cup 0.025 ml of the solution as the test compound 
was added by sterilised pipettes. The petridishes were kept in a cold room to facilitate 
the diffusion of the solvent for about 2 h. The plates were then incubated at 37°C for 
24 h. The extent inhibition was measured by the width of the zone in mm. Norfloxacin 
was used as a positive control and solvent control was also used to know the activity 
of the solvent (DMF). The results of the antibacterial testings have been summarised 
in Table 3. 


Table 3- Antibacterial and antifungal activity of biscoumannyl ethers 


Compound 

Zone of inhibition in mm 

Bacteria 

Psedomonas clororaphis 

after 24 h 

Fungi Aspergillus niger 

Ille 

18(++) 

22 (+ + +) 

Illf 

18(+ + ) 

14(+) 

IVa 

14(+) 

14(+) 

IVb 

14(+) 

14(+) 

Vd 

18(++) 

18 (+ +) 

Ve 

14(+) 

14{+) 


Antifungal Activity ; Fungicidal activity of test compounds were evaluated against 
Aspergillus niger species by cup plate method'®. The fungi were sub-cultured in the 
medium containing peptone (1%), yeast extract (0.6%), sodium chloride (0.5%), 
potassium dihydrogen phosphate (0.3%) and glucose (1%) in distilled water. The pH 
of the medium was adjusted to 6.0 and sterilised at 15 lbs for 20 min. The method of 
testing for antifungal activity was the' same as followed for evaluating antibacterial 
activity. Griseofiilvin was used as a positive control and solvent control (DMF) was 
also used to know its activity. The results of fungicidal activities are summarised in 
Table 3. 

Amongst the various biscoumarinyl ethers (HI) having R = 5, 6-Benzo and R = 6- 
C1 substitution in coumarin enhanced the antibacterial activity against Pseodomonas 
chlororaphis. Similarly in compound (V) having R=7,8-benzo substitution enhanced 
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the antibacterial activity. These compounds also reflected comparative anti-fungal 
potency. 
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Abstract 

The conductometric and electronic absorption spectroscopic studies of L-amino acids with 
iodine in 50% aqueous ethanol medium have been reported It is shown that the amino acids 
exhibit 1 :2 (ammo acid iodine) stoichiometry while 1 1 in excess of donor The band around 
212 nm has been assigned to T moiety while the bands near 293 and 355 nm are attributed to I's 
moiety The interaction of I's and 1'^ to donors is discussed 

(Keywords molecular interaction/charge transfer complexes/amino acids/iodme) 


Introduction 

Electronic absorption spectra of solutions of proteins namely fibrinogen, y-globin, 
trypsin, bovine plasma albumin, amino acids viz, glycine, tryptophan, hydroxyproline, 
tyrosine, a-alanine, (3-alanine, isoleucine, leucine, phenyalanine, arginine, and amines 
namely ethyl-, diethyl- and triethylamine with iodine which is essential constituent of 
thyroid hormone have been reported in the literature’. It has been demonstrated that 
the effect of amino acids and proteins in aqueous solutions of iodine can be 
reproduced by increasing the polarity of the solvent. The amino acids and the proteins 
exist as polar molecules in aqueous solution and the spectral changes observed on 
mixing in iodine has been attributed to the increase of the polarity of the solvent. 
However, aminobenzoic acids being non-polar molecules would not be expected to 
have such an effect unless the aminobenzoic acid-iodine interaction product is polar. 
No firm evidence has been produced that charge transfer occurs between the amino 
acids or proteins and iodine but in view of the results with the non-polar 
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aminobenzoic acids and chloranil' and the earlier evidence^"^ that the amino acids and 
proteins are good n-electron donors, it is likely that complexes of the form, (amino 
acids/or proteins, 1"^) are initially present in these solutions. Moreover, the equilibrium 
constant K which reflects on the stability of the complex and the enthalpy which is a 
true measure of the strength of interaction, however, has not been determined with any 
of the electron donors so far studied. 

We report in this paper the conductomertic and electronic absorption studies on 
the interaction of iodine with ten L-amino acid donors namely glycine, leucine, 
isoleucine, phenylalanine, tryptophan, valine, asparagine, threonine, histidine 
monohydrochloride and lysinemonohydrochloride. These studies are dimed at (i) 
establishing the stoichiometry and nature of CT complexes; (ii) determining the 
stability constants and thermodynamic parameters of the complexes; (iii) correlating 
the parameters characterizing the CT complexes with those of D and A, and (iv) 
substantiating or contradicting the earlier studies. 

Materials and Method 

The conductometric^’’ (Table 1) and thermodynamic parameters* (Table 2) were 
determined by experimental procedures and methodology reported earlier’. For 
composition determination, the equimolar stock solution (5 x lO"^ M) of D and A were 
prepared in 50 per cent aqueous ethanol whereas for thermodynamic paramerters, the 
stock solution of amino acids with concentration 5 x 10’’ M and acceptor with 5x10’^ 
M except in case of tryptophan -iodine system were prepared in 50 per cent aqueous 
ethanol. In case of tryptophan-iodine system, the concentration of D and A were 5 x 
10’* M and 5 x 10’* M, respectively Electronic absorption spectra of solutions of L- 
amino acids with iodine in 50% aqueous ethanol were measured on Backmann DU7 
Spectrophotometer. 


Results and Discussion 

The presence of a well developed maxima at a stoichiometric ratio of 1 :2 (amino 
acid : iodine) (Fig. 1) in the conductivity concentration plots observed in all the 
systems studied are indicative of the fact that in the various amino acids-iodine 
systems, molecular complexes of ionic CT type are formed. Inspection of data in 
Table 1 reveal that there exists no definite relationship between the values of Op, CTm 
and jE* among themselves as well as with pKi of amino acids as expected due to small 
variation of pKj (+0.34) and pK 2 (±0.74) values of amino acids”. 
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Table 1- Conductance (a) data for CT complexes of L-ammo acid donors with iodine acceptor in 50% 
aqueous ethanol at 293 K 


Donor 

pKi 

Stoichiometry 

Conductance (jiS) 

£:*'(kJ 


(amino 
acid) (a) 



0-0 

CTm 

mor') 

Isoleucme 

6 04 

(1.2) 

52 9 

41 3 

8025 

13 28 

Leucine 

6 03 

(1 2) 

62 5 

46 5 

1725 

8.85 

Phenylalanine 

5 67 

(1 2) 

52 6 

42 4 

8019 

7.29 

Tryptophan 

5 93 

(1 2) 

57 5 

45 5 

7535 

6 16 

Valine 

6,01 

(1:2) 

55.5 

44.0 

7468 

6 64 

Asparagine 

10.$4 

(1:2) 

58.5 

42 5 

10756 

8 45 

Glycine 

6 06 

(1-2) 

55 0 

44 6 

6662 

9.96 

Threonine 

5.59 

(12) 

57.5 

43 5 

9726 

8.12 

Histidmemonohydrochlonde 

7 76 

(12) 

59.1 

45.0 

8952 

9 96 

Lysinemonohydrochloride 

9 70 

(12) 

53 6 

43.0 

8217 

9 30 


60.01 


rAHO 
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The equilibrium constants of the complexes recorded in Table 2 have been 
determined on the assumption that in large excess of the donor ), the 

composition of the complex is 1: 1. The correctness of the assumption has been 
verified fi'om the consistency of the results obtained from model^ based on 1:1 
complex formation. The absence of any definite relationship between the parameters 
of Tables 2 i.e. K and pd with those of Table 1 i.e. Op, ^ and pKi may also be due to 
the reason cited earlier. The order of stability of the complexes (i.e. K value) is: 
histidinemonohydrochloride > lysinemonohydrochloride > threonine > glycine > 
phenylalanine > tryptophan > isoleucine > valine > leucine > asparagine. On the other 
hand the enthalpy (-A/7) which i^ a measure of the molecular interaction follows the 
trend: tryptophan > leucine « phenylalanine > theronine > isoleucine > glycine > 
asparagine valine > histidinemonohydrochloride > lysinemonohydrochloride. The low 
stability {K) and strong interaction (-A//) of tryptophan as compared to that of 
phenylalanine is because of imidazole N in the former. 

Further, the inspection of Table 2 (columns 8-10) shows that pa the depolarization 
constant grows with the growth of temperature which is reverse to the behaviour of K 
values with temperature as expected*. 

The absorption spectra (Table 3) of dark brown colored old solution of iodine in 
de-ionized water exhibits bands at 287, 353, and 465 nm, the first two bands being 

characterized of the tri-iodide negative ion (Ij) formed according to the following 
reactions”. 

I2 r + r 

h + r^==^ I3 

The 465 nm band arises from solvated iodine. The iodine absorption band in inert 
solvent is 540 nm'^. The absorption spectra of iodine in potassium iodide solution 
(freshly prepared) in water exhibit bands at 230, 285.5 and 352.5 nm (Table 3). The 

later two bands may be attributed to the tri-iodine ion ( I 3 ) as expected. The band at 
230 nm may be because of iodide (T) as expected”. When iodine is directly dissolved 
in neat alcohols the conductance increases with time and reaches the maximum value 
in about 24 h”. It is reported that the intensity of CT band at 225 nm due to T in 
iodine-ethanol system decreases with time and also progressive intensification of I's 
ion absorption at 297 and 363 nm have been observed”. Thus, it can be inferred that 
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time dependence of electrical conductivity is due to the transformation of initially 
formed 1 : 1 outer charge transfer complexes followed by veiy fast reaction of the inner 
charge transfer complexes with iodine to form the tri-iodide ion 


Tables- Electronic absorption spectra of L-amino acids-iodine systems in 50% aqueous ethanol 
solution 


System 

^max (Absorbance) 

Iodine m deionised water 

287, 353, 465 

Iodine in 50% aqueous ethanol 

213, 226, 290 5,259 5 

Kl+Iodine in 50% aqueous ethanol 

230, 285 5, 352 5 

L-Amino acids+ Iodine (I 2 ) in 50% aqueous ethanol 


Isoleucine 

212 5 (0 5672), 297 (.2482), 356.0 (0 2498) 

Leucine 

212 0( 6301), 297 0( 3699), 359.0 ( 3695) 

Phenylalanine 

212 0( 4957), 294 0( 2338), 357.0 ( 2325) 

Tryptophan 

212 5(1 0345), 293 0 ( 0835), 347 0 ( 8040) 

Valine 

212 0 ( 9917), 291 0 ( 1760), 358 5 ( 101 1) 

Asparagine 

212 5(1 0192), 295 5 ( 2514), 358 5 ( 1387) 

Glycine 

212.5( 7072), 293 0(.1698), 359 0 ( 1037) 

Threonine 

212.5 (1 0981), 291 0 ( 3345), 358 5 ( 1918) 

Histidinemonohydrochloride 

212 0 (2 0152), 288 5(0 2425), 359 0 (0 1308) 

Lysinemonohydrochlonde 

212.0 (2 0131), 299 5 ( 4181), 356.5 (.4191) 


very fast 

ROH + h ^ ROH - 12 ► ROH Tl" 

Outer charge Inner charge 

transfer complex transfer complex 

fast 

ROH rr+i 2 X ROH rri 2 ►roh r+ij 


R = CH3, C2H5,n - C3H7, i - C4H8 and CeHsCHj 
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In case of amino acids-iodine systems in 50% aqueous ethanol (Table 3) the last 
two peaks (291-299, 347-354 nm) are assigned to Ij whereas the peak around 212 nm 
may be attributed to T since the cc nduct^' metric studies give rise to a complex having 
stoichiometiy 1:2 (D:A). The tej- ative Jtructure of the molecular complex may be 
thus written as: 

o 

0 

/ 

nh^-ch2-c^ r 



The interaction of Ij to protonated N in the present studies finds support from 
earlier spectroscopic observations that the band at 226 nm in the absorption spectra of 
iodine with amines, amino acids and protein in aqueous solution is due to binding off 
with protonated nitrogen. 
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Abstract 

The stability constants of Zn^*, Co^\ Ni^* and Mn^* complexes of 3-(4- 

methoxyphenyl)-l-(2’-hydroxynaphthyl)-2-propen-l-one (HMPHPO) have been determined in 
75% (v/v) dioxane-water, acetone-water and 2-ethoxyethanol-water mixtures using the 
Bjerrum-Calvin pH titration technique as developed by Irving and Rossotti. The stability 
constants are in fair agreement with Irvmg-Williams order The results obtained are compared 
with the data available in the literature Abnormal behaviour of the chelates was observed in 
acetone-water medium. The thermodynamic parameters such as free energy, enthalpy and 
entropy change involved in complexation have been calculated 

(Keywords 2-hydroxy- 1-acetonaphthone / p-methoxybenzaldehyde/free energy/stability 
constants/ enthalpy/entropy/Calvm-Bjerrum technique) 

Introduction 

Very few systems have been reported' showing the relationship between thermal 
stability of metal chelates and structure of the chelating agent. Wendlandt and co- 
workers^'^ and Hill and co-workers ’ studied the thermal properties of metal chelates 
with different types of complexing ligands. Lense et al.^ reported that o- 
hydroxychalcones are much more reactive with metal ions than the ketones and 
aldehydes from which they are synthesized. Holm and O’Connor® have reviewed 
structural studies on several metal chelates of P-diketones and 2-hydroxycarbonyl 
compounds. The metal chelates of divalent transition metals with some o- 
hydroxychalcones standout as a distinct class of o-hydroxycarbonyl compounds with 
low spin square-planar configuration, which do not easily form adducts. Extensive 
conjugation was found to he responsible for the strong field nature of the ligand'®. 
Chalcones were found to be physiologically and pharmacologically active. They were 
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also found to have germicidal", bactericidal'^ fungicidal’^ and carcinogenic 
activity". Several of the chalcones have been employed as useful analytical reagents. 

The present investigation deals with the determination of stability constants and 
the thermodynamic parameters of the complexes of UO 2 , Cu , Zn , Co , Ni and 
Mn^"' with 3-(4-methoxyphenyl)-l-(2’-hydroxynaphthyl)-2-propen-l-one(HMPHPO, 
Fig. 1) at different temperatures employing Calvin-Bjerrum pH titration technique"'" 
as modified by Irving and Rossotti’®. 



Fig l-StTUcture of HMPHPO 


OCHs 


Materials and Method 

The ligand 3-(4-methoxyphenyl)-l-(2’-hydroxynaphthyl)-2-propen-l-one 
(HMPHPO) was prepared by the condensation of 2-hydroxy- 1-acetopaphthone with p- 
methoxy benzaldehyde employing the Claisen-Schmidt condensation reported 
earlier". The purity of the sample was checked by elemental analysis and melting 
point determination. The solution of the ligand was prepared using dioxane, acetone 
and 2-ethoxyethanol as solvent. The metal nitrate solutions were prepared in CO 2 free 
doubly distilled water and were standardized by standard methods. Dioxane, acetone 
and 2-ethoxyethanol were purified before use. An ELICO (< 1-120) digital pH meter 
was used for pH measurements. 

The experimental method consists of pH titration of the following carbonate 
free solutions of 50 ml at the desired temperatures against standard sodium 
hydroxide. 
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(a) 5 ml 0.1 M KNO3 + 37.5 ml organic solvent + 5 ml 0.01 M HNO3 + 2.5 ml 
water. 

(b) 5 ml of 0.1 M KNO3 + 1 ml of 0.1 M ligand in pure organic solvent + 36.5 
ml organic solvent + 5 ml 0.01 M HNO3 + 2.5 ml water. 

(c) 5 ml of 0.1 M KNO3 + 1 ml of 0.1 M ligand solution + 36.5 ml organic 
solvent + 5 ml 0.01 M HNO3 + 0.5 ml of (0.01 M) metal solution + 2 ml 
water. 

The pH values (B) were corrected in all the aqua-organic mixture by the method 
of Van Uitert and Hass^° (- log [H^ = B + log Uh) and others^' under the present 
condition the following logUn values were obtained in different media under 
consideration : 0.36 (75% aq. Dioxane), 0.29 (75% aq. acetone) and 0.24 (75% aq. 2- 
ethoxyethanol). 


Results and Discussion 

In our earlier work^^, calculation of proton ligand formation constants of the 
ligand was explained, the same method is applied for this ligand (HMPHPO) and 
calculated proton ligand formation constant in 75% aq. dioxane at 30®C, 35®C and 
40'’C are 12.80, 12.54, and 12.41 respectively. In 75% aq. acetone and 75% aq. 2- 
ethoxyethanol at 30°C are 12.36, 11.77 respectively. 

From the titration curve of solution of (b) and (c), n and Pl values were 
calculated’^’ The stepwise formation constants were corrected by least squares 
method'^ and are presented in Table 1. 

Table 1- Stability constants of bivalent transition metal complexes with HMPHPO at different 
temperatures and different solvent media (Ionic strength 0 IM) 





75% Dioxane 

75% 

Acetone 

75% 2- 
Ethoxyethanol 

sox 

35X 

40 X 

SOX 

30X 


\Ogki 

10.98 

10.74 

10 56 

10 73 

10.38 


log *2 

10.36 

10.17 

9.78 

10 34 

10 10 


log Pz 

21 34 

20 91 

20 34 

21.07 

20.48 


Table 1 Contd... 
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Table 1 Contd... 



log*! 

9 86 

9 64 

9 48 

961 

9.52 

Cu^-^ 

log *2 

9 43 

9 32 

9 15 

9 35 

8 85 


log ^2 

19.29 

18.96 

18 63 

18 96 

18 37 


log A:, 

9 08 

8 83 

8 55 

8 74 

8 55 

Zr?-* 

logkj 

8 36 

8 12 

7 94 

8 42 

7 96 


log P 2 

17 44 

16 95 

16 49 

17 16 

165 


log ^1 

8 35 

8 01 

7 88 

8 02 

7 82 

Cc?* 

logi2 

7 64 

7 58 

7 45 

7 48 

701 


log P 2 

15 99 

15 59 

15 33 

15 50 

14 83 


log fei 

8 05 

7 78 

7 55 

7 94 

7 72 


logA:2 

7 48 

7.35 

7 24 

7 38 

6 84 


logP2 

15 53 

15 13 

14 79 

15 32 

14 56 


log A, 

7 82 

7 66 

7 48 

7 62 

7 34 


iogk2 

7 04 

6 82 

6 64 

7 08 

661 


log P 2 

14 86 

14 48 

14 12 

14 70 

13 95 


The formation of metal chelates takes place by replacement of hydrogen of the 2- 
hydroxynaphthyl groups and coordination occurs through the oxygen of the ketooic 
group. The acid dissociation constant of the ligand is found to decrease m different 
media m the order dioxane > acetone > 2-ethoxyethanoL A plot of pkL against HD is 
linear. However in acetone-water solvent the dissociation constant of the ligand is less 
than the expected value and it may be attributed to the non-ideal behaviour of acetone- 
water solvent. The order of stability of the metal chelates, observed and studied in all 
the three mixed solvents is > Cu^"^ > Zn”"^ > > Ni^"^ > Mn^"^ This order of 

stability follows the Irving - Williams^^ natural order expect cobalt (II) leading to the 
conclusion that the ligand used in the present investigation is weak field one. The 
greater stability of cobalt (11) chelate when compared to nickel (II) chelate may be 
attributed to the additional stabilization due to John - Teller distortion present m the 
former case. The order cobalt (H) > nickel (II) has been reported for some other 
ligands also viz. 4-nitro-2-aminophenoI^'^, 8-hydroxy-2-methylquinoIinol^^, 2- 
hydroxychalcone^^ and 3 -(4-chIorophenyl)- 1 -(2 ’-hydroxynaphthyI)-2-propen- 1 -one^^. 
The values of log P for zinc (II) chelates are greater than those of nickel (II). The 
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order zinc (II) > nickel (II) has been observed earlier also^’’ where the relative low 
value for nickel (H) has been attributed to the presence of a square planar structure. 
Further the hydrolysis constant of zinc (II) chelate is greater than that for nickel (II) 
indicates a greater affinity of zinc (II) for - OH ions. This factor seems to be 
accurate for the order zinc (II) > nickel (II). The order of stability of complexes with 
respect to the solvent in dioxane > acetone > 2-ethoxyethanol. The dielectric constants 
for the pure as well as mixed solvent are in the order acetone > 2-ethoxyethanol > 
dioxane. This also shows the abnormal behaviour of acetone-water mixture, similar 
abnormality was reported in the literature^^’^^. A plot of log p vs. ionization potential 
shows that the stability of the complexes can be better correlated with the second 
ionization potential rather than the first. 

The data presented in Table 1 shows that the difference in log ki and log ki values 
are small and the ratio of log ki / log ki is positive in all cases indicating that there is 
little (or) no steric hindrance to the addition of the second ligand molecule^^. It is also 
evident from the values of log k\ and log ^2 that there is a almost equal tendency for 
the formation of neutral complex species ML 2 as for the species Ml"^. 

The stability constant values decrease with increase in temperature (Table 1) that 
is in agreement with the conclusion reached by Pitzer^*. The thermodynamic 
parameters free energy (AG), enthalpy (AH) and entropy (AiS) change have been 
calculated^® from the stability constants of the metal complexes determined at three 
different temperature (30° C, 35° C and 40° C) in 75% dioxane and are given in Table 
2 . 


Table 2- Thermodynamic parameters of transition metal complexes with HMPHPO in 75% dioxane 


Metal 

-AG 

(K J.mole-') 

-AH 

(K J mole') 

4-A5 

(J deg'^mole*^) 

uo^ 

123 79 

113 02 

35.54 


111 90 

107.62 

14 13 

Zx?* 

101 17 

94 19 

23.12 

Co^^ 

92 76 

91 88 

2.11 

Ni^ 

90 09 

80 71 

30 96 


86 20 

70.66 

51 28 


AG values at 30°C 

Error limits are ± 0.01 for AG and AH and ± 0.1 for hS 
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The AG values are all-negative and suggest that complexation is a spontaneous 
process. These values become less negative at higher temperatures indicating that low 
temperature favours complexation. The negative A/f values for complex formation 
indicate the release of energy during the interaction of metal ions with ligand. The A/T 
values of the metal complexes (HMPHPO with different metal ions) follow the order 
> Cu^"^ > Zn^^ > Co^"^ > > Mn^^, in agreement with the natural order of the 

stability constants. The values of AS are positive for all the system indicating that 
entropy is favourable for the complex formation. 
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Abstract 

QSAR studies of anti inflamniatory(l-lO) novel tetrahydro 2H-isoindoles have been discussed 
The inhibitory potency of these compound is found to be dominantly controlled by steric factor 
'x* Significant correlation have been obtained between inhibitory potency and molecular 
connectivity 'x*. 

(Keywords QSAR/1,3 diaryl-4,5,6, 7 -2 H-isoindole/Cox-2) 

IntroductioB 

Kujubu et al? suggested that an inflammatory or mitogenic stimulus might result 
in gene expression responsible for the synthesis of a second inducible cyclo 
oxygenase isoenzyme called Cox-1 Cox-2^ rapidly emerged as an invaluable target for 
developing anti inflammatory agents with reduced gastrointestinal or renal secondary 
effect .Nonsteroidal anti inflammatory drugs (NSAIDS)are widely used to treat pain, 
fever and inflammatory condition including Osteoarthritis. Recently ,it has been 
postulated that Cox-2 inhibitor may effectively relieve symptoms in diseases such as 
Osteoarthritis, while exhibiting a considerably safer toxicity profile than classical 
NSAEDS especially with respect to gastric complications^. It has been recently 
demonstrated that IL-l-stimulated chondrocytes from normal cartilage expressed the 
cox-2 gene and produced the cox-2 protein followed by PGE2 release'*, Cox-2 protein 
over expression has been detected in chondrocytes from Osteoarthritis patients as 
compared with normal subjects whereas cox-1 was not detected in either normal or 
Osteoarthritic chondrocytes^. PGE2 over induced in Osteoarthritic cartilage has been 
shown to coincide with cox-2 impregnation hence suggesting that PGE2 may be 
differentially regulated in normal and Osteoarthritic cartilage*.In this field only 
Nimesulide’, its analogues NS 398* and structural flosulide analogues®’*® are described 
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as preferential cox-2 inhibition". Pyrazole derivatives by Searle are serious candidates 
for potent and selective inhibitors which were structurally related to the class of 
cyclopentenyl lactones*^’^^ Celeocoxib and rofecoxib(vioxx, 2) have been 

recently marketed. 



After several modification in structure in order to enhance the cox-2 inhibition 
potency a lead and a compound (3) had been reported 



Method 


In present work steric parameter i.e. molecular connectivity widely has been 
used to assess the role of molecular size, with a view to provide a better rational 
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approach for the design of potent drugs .The (first order of molecular connectivity 

index) was calculated as suggested by Kupchiik**''®. 

The regression analysis was performed with the help of statistical package SPSS. 

The data listed in Table 1 has been used for QSAR study. ICso, Cox-2 refer to 
concentration of compounds required to achieve 50% inhibition against Cox-2 from 
mouse resident macrophages 


Table I- Physio-chemical parameters and biological activity of Tetrahydro-2H isoindole. 


s. 

No. 

Structure 


Physio-chemical 

parameters 

Activity-log I/IC 50 


R 

Ri 

Rz 

X&q 

Id 


Obs 

cal eqn 
( 1 ) 

cal eqn 
( 6 ) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 . 

H 

H 

H 

2.313 

4.418 

7 528 

8.824 

8.798 

8.673 

2 . 

H 

H 

H 

2.313 

4417 

7.428 

8.481 

8.870 

8.743 

3 

NH2 

H 

H 

2.308 

6.450 

7.984 

8 744 

8 469 

8.357 

4. 

NHSO2CH3 

H 

H 

2.339 

4.685 

10.786 

6.301 

6 444 

6.412 

5. 

H 

F 

F 

2.363 

4.518 

7 978 

8 869 

8.473 

8.798 

6. 

H 

CH3S 

CH3S 

2.331 

5.013 

11.162 

6.301 

6 173 

6.152 

7. 

H 

CH3 

CH3 

2.290 

4518 

8.599 

7 778 

8.028 

7.930 

8 . 

H 

OCH3 

OCH3 

2 323 

5 013 

8.825 

7.671 

7 862 

7.774 

9 . 

H 

Cl 

Cl 

2.339 

4.518 

8.734 

8 301 

7.927 

8.273 

10 . 

H 

F 

C3H3N2 

2 353 

4.757 

9.314 

7.377 

7.508 

7.435 


The regression analysis of the data presented in Table 1 gave following significant 
equation: 

Log(l/ICjo) = -0.722(±0 07) 'x* + 19.23 (±0.64) (1) 

« = 10, r = 0.97, =0.93, J= 0.25, F=100.9 
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where n represents the number of data presents, r the correlation coefficient, F the 
varience ratio between the calculated and observed activities and s the standard error 
of estimate. Equation (1) gave significant but negative coefficient accounts for 92.7% 
of varience and F ratio was found to be significant at 99% level [(i^ 2 ,io) = 7.8] 

In order to study the effect of specific substituents at R, R\, Rj position six 
indicator parameters were taken, Ind-i for compound possessing NH 2 at pyrrole ring, 
Ind -2 for compound possessing NHSO 2 CH 3 at position R, Ind -3 for compound having 
F at Ri and R 2 position, Ind .4 for compound CH 3 S,Ind _5 for compound OCH 3 and Ind _6 
for compound Cl and F at position R\ and R 2 phenyl ring. Following significant 
equations were obtained. 

Log(l/IC5o)=-0.703(± 0.07) 'x* +0.324(±0.294)Indi + 14.03 (2) 

« = 10, r = 0.968, P = 0.937, 5 = 0 27, F = 52.39 

Log(I/IC5o)=-0.692(±0.08) ‘x* -0.224(±0.35)Ind2+ 13.99(±0.77) (3) 


n = 10, r =0.97, =0.93, s = 0.25, F=46.8 


Log(l/IC 5 o) =-0.702 (±0.07) ’x* +0.348(±0 29)Ind3+ 14.02(±0.64) (4) 


« =10, /■ =0.97, r^=0.93, s=0.26, F = 55.93 


Log (1/IC5o)=-0.761 (±0.09) ‘x* +0.24(±0.39)Ind4 ± 14.54(±.84) (5) 


« =10, r = 0.964, r^=0.930, s=0.28, F = 46.6 

Log(l/IC 5 o) = 0.722(±0.07) ‘x* - 0.21 1 (±0.29)Ind5 +14.25(±0.66) (6) 

n=l0,r =0.97, =0.93, 5=0.28, F = 47.57 

Log(l/IC5o)=-0.694(±0.05) ’x* +0.436 (±0.17)Ind6± 13.89(±0.51) (7) 


n =10, r =0.90, r^=0.97, s=0.27, F =55.42 



quantitative structure activity relationship studies 


119 


Equation(2-7) account for 93.9%, 93.0%, 93.1 % variance of F are also found to be 
significant for all above equations at 99% level [(/='2,io)=7.56]. 



Fig 1 



Fig 2 







120 


ARUN K. SRIVASTAVA et al 


Conclusions 

On the basis of the above study it can be inferred that steric factor with 
negative (-) coefficient plays an important role in determining the Cox -2 inhibitory 
potencies-the negative coefficient shows that with small increase in constituent size 
potency decreases sharply. The coefficient value of Indi, Inds and Inde are having 
positive (+) coefficient suggesting that -NH2 at pyrrole ring (Indi), Halogen (Cl or F, 
Inds and Ind6) particularly F at para position of phenyl ring shows higher inhibitory 
potencies. Whereas bulky substituent like -NHSO2SO3 at pyrrole ring Ind2 and -OCH3 
at para position (Inds) shows decrease in inhibitory potencies however substituent 
~CH3S Ind4 shows positive coefficient and increase inhibitory potencies suggest there 
may be an active pocket for CH3S where S (Sulphur) seems to be playing active role, 
therefore smaller substituents group like halogen, amine at Ri and R2 position at 
phenyl ring should be preferred in future drug synthesis. In the present series these 
findings are in conformity with Portevin et. al\ These fact are also proved by Fig. (1 
& 2) plotted between predicted and calculated inhibitory potency, steric parameter and 
log I/IC50 respectively. 
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Abstract 

The paper deals with the problem of the reflection and refraction of plane waves at an interface 
between two dissimilar micropolar elastic solid half-spaces with stretch The amplitude ratios 
for various reflected and refracted waves are calculated theoretically Numerical values of the 
amplitude ratios are computed and plotted against the angle of emergence of incident 
longitudinal microstretch wave 

(Keywords reflection/refraction/micropolar elastic solid/axial stretch/amplitude ratios/ LMS 
wave) 


Introduction 

It is well known that material response to external stimuli depends heavily on the 
motion of its inner structures. Classical elasticity ignores this effect by ascribing only 
translation degrees of freedom of material point of body. Micropolar theory of 
Eringen and his coworkers’’^ by including intrinsic rotations of the microstructure, 
provides a model that can support body and surface couples and display high 
frequency optical branch of the wave spectrum. 

Eringen^ extended his work to include the effect of axial stretch during the 
rotation of molecules and developed the theory of micropolar elastic solid with 
stretch. For engineering applications, composites with rigid chopped fibres, elastic 
solid with rigid granular inclusions, and other industrial materials such as liquid 
crystals may be modeled as micropolar elastic solids with axial stretch. 
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During last four decades, various researchers'*'''* discussed some problems on 
waves and vibrations in micropolar elastic solid. Recently, Tomar and Gogna'®" 
discussed some problems on reflection and refraction at an interface between two 
micropolar media. The present investigation is concerned with the reflection and 
refraction of plane waves at an interface between two dissimilar micropolar elastic 
solids with stretch. 


Basic Equations and their Solutions 

Following Eringen^’^, the constitutive and field equations in a micropolar elastic 
solid with stretch without body forces and body couples can be written as 



(1) 

mkl=^0£rk\^\r + a (|)r,r5il+p(|)*,l+Y(t)l.*. 

(2) 

pi= ao<I>*,t+(po/3)£rM<|)M, 

(3) 


and 


(cf +C3 )V (V. u)-(c2 +c|) )V X (V X m)+ cl'Vx^^ii ( 4 ) 

(C4 +C5))V (V. (j))-C4 V X (V X (j)) + (olV X u - 2 a)f <|) = <i), ( 5 ) 


C6^V^<D*-riO*=(J)*, 

where 

cf=(X.+2p)/p, 

cl=y/pj, 

cl=lao/pj 


C2=p/p, 

c^=(a+p)/p7, 


ri=2Tio/ft/. 


( 6 ) 


C3=K/p, 

=Cj/j=K/pj, 

(7) 



WAVE PROPAGATION IN A MICROPOLAR ELASTIC SOLID 


125 


where symbols X, k, a, P, y, ao, po, t|o, p have their-usual meanings, u, (j) and O* are 
displacement vector, microrotation vector and scalar micro stretch respectively. 6 */ is 
the kronecker delta. The superposed dots on right hand side of equations (4) to ( 6 ) 
denote the second order partial derivatives with respect time. 

Parfitt and Ermgen^ have shown that there exists four basic waves propagating in 
an infinite micropolar solid, namely, a longitudinal displacement wave (LD wave) 
travelling with speed Vi, a longitudinal microrotational wave (LM wave) travelling 
with speed V 2 and two sets of two coupled transverse displacement and 
microrotational waves (CD I and CD II waves) with velocities Vs and F 4 . 

In order to solve equation ( 6 ), we assume that wave propagates in a positive 
direction of unit vector n as 


O* = D exp[i^(n.r-F/)], ( 8 ) 

where D is the amplitude of the wave, V is the phase velocity, r is the position vector, 
k is the wave number and co = kV. 

Substituting equation ( 8 ) in equation ( 6 ), we get 

cl )- - ri]D exp[ 2 A:(«.r - V()] = 0, (9) 

where ri = rio/oto- 

From equation (9), we have 

V^= cll{\-{_r2 cllis,'^)}, (10) 

Equation ( 10 ) gives the speed of the wave given by equation ( 6 ) and we call the wave 
as longitudinal micro stretch (LMS) wave. 

V will be real and finite only if 


l-('‘2 Cg/CO^) > 0, 
CO > (2 rio/py)’^^ , 


( 11 ) 
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This is the condition for the existence of the longitudinal micro stretch wave. 

Thus, in an unbounded micropolar elastic solid with stretch there exists five basic 
waves travelling with distinct speeds. 

Since we are discussing a two dimensional problem in x-z plane, we have 


U =(MiO,«3), 


where 


3 ^ - 5 ^/ 
dx dz 


0 ), 


_ 3^ 3v|; 

^ dz dx 


( 12 ) 


(13) 


where <j), vi/; are potentials satisfying certain wave equations in micropolar media with 
stretch. 


Boundaty Conditions 

For two-dimensional motion in xz-plane, the appropriate boundary conditions at 
the interface z = 0 are the continuity of normal stress, tangential stress, couple stress, 
normal vector first moment, tangential vector first moment, normal displacement, 
tangential displacement and microrotation, i.e., 


II 

N 

/ = / ' 

*ZJC 5 

II 

s 


P'., 


= P/ 

W3 = U3\ 

Ui = Ui\ 

<1>2 = 

^ 2 , at z = 0, 

(14) 


The quantities with primes correspond to the medium (z < 0) and those without primes 
to the medium (z > 0) throughout this manuscript. 

Reflection and Refraction 

We consider two micropolar elastic solid half-spaces with stretch which are in 
welded contact along a plane interface z = 0. We introduce rectangular Cartesian co- 
ordinates (x, y, z) and place the origin at the interface separating the two half-spaces. 
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We now consider a plane harmonic seismic body wave [LD, CD I or LMS] with time 
dependence proportional to exp(/coO, propagating through the medium (z > 0) and 
striking obliquely at the interface z = 0 and making an angle 0o with the interface, 
Corresponding to this incident wave, we get waves in medium (z > 0 ) as reflected LD, 
CD I, CD II and LMS waves and refracted LD, CD I, CD 11 and LMS waves in 
medium (z < 0). 

The non-vanishing potentials suitable for the problem are given by 

^ = Aq exp[/^i(x cos 6q-z sin 6o)+/o)i/] 

+A] exp[/A:i(jc cos 0 i+z sin 0 i)+/<Bit], ( 15 ) 

vj/= Ao exp[/ 5 i(x cos do-z sin 0o)+/co2t] 

+A2 exp[z8i(x cos 02+r sin 02)+/co2^]> 

+A3 exp[/52(x cos 03+z sin 03)+/©3t], ( 16 ) 

([>2= EAo exp[i8i(x cos do-z sin 0o)+/a)2^] 

+£A2 exp[/ 5 i(x cos 02+z sin 02)+z©2?] 


+FA3 exp[z82(x cos 03+z sin 03)+/(fl3t] ( 17 ) 

<!>*= GAo exp[/i4(x cos Oq-z sin 0o)+/®4t] 

+GA4 explii:4(x cos 04+z sin 04)+/co4t], ( 18 ) 

([)'= A(i'exp[/A:i'(x cos 0 i'-z sin 0 i')+z©i't] ( 19 ) 

\|/'= -42'exp[i^i'(x cos 02'-z sin 02')+i®2'O 

+A3'exp[iS2'(x cos 03'-z sin 03')+/©37], (20) 
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<t) 2 '= E'Ai'txpliSi' (x cos 02-2 sin 62 0 +*fi> 2 '^] 
+F 'A 3 'exp[/ 82 ' (x cos 03 '-z sin 03 't\, 

<!>*= G'A^'txplik^' (x cos 04 '-2 sin 04 ')+iO 4 '/] 

where 

E= 8f(8^ -^-pq)/p(2q~^) 

C2 +C 3 C 4 

F= 8^(8^ — 2^-T-P9)/M2g-'^) 

C 2 +C 3 C 4 

where 


( 21 ) 

( 22 ) 


(23) 


p = ii/(p+it), q=K/y, 8 ? = X>^ 5l= X\(i)\ ± (5MaC)'^^]/2, 


5 = 


g(p- 2 ) , 1 1 

CO (C2 + Cj ) C4 


C = 


1 


(C 2 +C 3 ) 


1 2 q 


^-4 


CO 


(24) 


and G is the constant of dimension L'^. Similar expressions with primes correspond to 
medium (z<0). Ao_ Ai, A 2 , A 3 , A 4 , Ai', A 2 , A 3 ' and A 4 ' are the amplitudes of incident, 
reflected and refracted waves. For incident LD wave, yfo=0 in equations (16) to (18); 
for incident CD I wave, yfo 0 in equations (15) and (18) whereas for incident LMS 
wave, Ao = 0 in equations (15) to (17). 

The extension of Snell's law will be 


F*'' cos 0 = Ff* cos 01 = A ,1 cos 02 = A, 2 'C 0 S 03 = F’ cos 04 


= Fj'"'ops'0i = Xi' cos 02 = X 2 ' cos 03 = F’ cos 04 


(25) 



W^VE PROPAGATION IN A MICROPOLAR ELASTIC SOLID 


129 


and 

coi=(02=co3=a)4=coi'=co2'=<03'=<B4'=<o(say) 


(26) 


where F* = Fi, F 3 or V, for incident LD, CD 1 or LMS wave respectively. 

Making use of the potentials given by equations (15) to (22) in the boundary 
conditions (14), after using the equations (1) to (3), (12) and (13), we obtain a system 
of eight nonhomogeneous equations as 

'Za,jZj=b„ (i=\, 2 , ,8) (27) 

7=1 


where Zj are the amplitude ratios for reflected and refracted LD, CD I, CD n and LMS 
waves, a,j are the coefficients of amplitude ratios and the relation between b, and a,j 
are given as follows 

(a) for an incident LD wave 


bi — ci\\, bz — azi, bi — 031, 64 - 041, bs — Osi, 


be - ^61) bz <371, bz - <3gi. 

(b) for an incident CD I wave 


^1 ~ On, bz — azr, bz - a^z, ^4 — ^i:42> be — a^z. 


be — 0625 bz — 072 , b% — og2 

(c) for an incident LMS wave 

b\ = — oi4, bz — O24, bz = — O34, 64 = 044, be — — O54, 
be ~ 064, bz = — O74, bz = og4 

If we neglect the stretch effect, we can obtain corresponding results for micropolar/ 
micropolar interface as shown by Tomar and Gogna*'*’” 
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Numerical Analysis 

Here, we shall confine our study for the case of incident LMS wave only. 
Following Gauthier'^, we take the following values of relevant parameters for 
micropolar elastic solid with stretch ( medium z >0 ) as 

p = 2.19 glcm\ X = 7.59 x lO” dyneW, 

p = 1.89 X lO" dyne/cm^, 

y= 0.0268 X 10“ dyne, j = 0.0196 cm^ k=0.0149 x lO" dyne/cm^ 

tto = 0.0915 dyne, Po = 0.0726 dyne, 

r|o = 0.0532 dyne, G - 0.05 cm’^. 

and, for medium (z < 0), we consider 

p' = 2.06 g/cm^, 

p' = 1 .62 X 1 0” dyne/cm^, 

y’ = 0.0254 X 10"^e, 

k' = 0.0134 X lO" dyne/cm^, 

ao' = 0.09 1 1 dyne, po' = 0.07 1 7 dyne, 

Tio' = 0.0528 dyne, G' = 0.045 cm'l 

For the above values of relevant physical constants, the system of equations (27) is 
solved for amplitude ratios by using Gauss-elimination method for different angle of 
emergence varying from 0° to 90° when coV cof = 10. The variations of these amplitude 
ratios with the angle of emergence of incident LMS wave are shown in Fig. (1 and 2). 


X' = 6.8 X lO” dyne/cm^. 


/ = 0.0187 cm^ 
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In Fig. 1, the variation of the amplitude ratios for reflected LD-wave, coupled 
waves (CD I and CD II) and LMS wave are shown by lines S,D(1), D(2) and D(3) 
respectively, whereas the lines S, D(l), D(2) and D(3) in Fig. 2. show the variation of 
the amplitude ratios for refracted LD-wave, coupled waves (CD I, CD II) and LMS 
wave respectively. 

The variation of the amplitude ratios for reflected and refracted LD waves are 
observed oscillatory in nature. The variation of the amplitude ratios for reflected 
coupled waves and refracted coupled waves are observed similar. They first increase 
to their respective maxima and then decrease. The amplitude ratio for reflected LMS 
wave, first decreases to its minima near Op = 89° and then increases sharply to its 
maximum value. The amplitude ratio for refracted LMS wave, first increases to its 
maxima near 0o = 89° and then decreases sharply to its minimum value zero. The 
amplitude ratio for reflected and refracted LMS waves are shown in Fig. (1 and 2) by 
D(3) after multiplying their original values by 10"*. 



Fig. 1-Variation of amplitude ratios for reflected waves with the angle of 
emergence of incident LMS wave. 
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Fig. 2-Vanation of amplitude ratios for refracted waves with the angle of 
emergence of incident LMS wave 

Conclusions 

The basic equations of motion for a micropolar elastic solid with stretch are 
solved. It is found that there exists a new wave, which we call longitudinal 
microstretch wave (LMS wave) in addition to waves shown by Parfitt and Eringen^. 
Theoretical calculations are presented for the cases of LD, CD I and LMS waves, 
incident at the interface between two dissimilar micropolar elastic solids with stretch. 
Numerical results are presented graphically for incident LMS wave only. If we neglect 
the stretch effect, the new wave (i.e LMS wave) will not exist. The problem though 
theoretical but is of physical interest in the fields of seismology, geophysics, 
earthquake engineering, etc. 
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Abstract 

The present paper deals with the study of hydromagnetic free convection and mass transfer 
flow of an incompressible, electrically conducting, viscous fluid past an infinite vertical porous 
plate embedded in porous medium of time dependent permeability under oscillatory suction 
velocity normal to the plate. It is considered that the uniform magnetic field is applied normal 
to the flow and the permeability of the medium fluctuates with time. The velocity distribution, 
temperature distribution, concentration distribution, fluctuating parts of the velocity, amplitude • 
and phase lead of skin-friction, rate of heat and mass transfer have also been obtained The 
behaviour of important parameters on the results obtained are discussed with the help of graphs 
and tables. 

(Keywords . free convection/mass transfer/porous medium/variable permeability). 

Introduction 

The phenomenon of free convection arises in the flow of fluid when temperature 
changes due to density variations leading to buoyancy forces acting on the fluid 
elements. Such flows have been the object of extensive research due to their enhanced 
concern in science and technology about buoyancy induced motions in the 
atmosphere, in bodies of water and quasi-solid bodies such as earth. Unsteady 
oscillatory free convective flows play an important role in chemical engineering, 
turbomachinery and aerospace technology. In addition, many transport process exists 
where the transfer of heat and mass takes place simultaneously as a result of combined 
buoyancy effects of thermal diffusion and diffusion of chemical species. 

Several authors*"® have studied free convection and mass transfer flow of a 
viscous fluid through porous medium assuming the permeability of the porous 
medium to be constant. However, the porosity of the medium may not necessarily be 
constant because the porous material containing the fluid is a non-homogeneous 
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medium. In the light of this fact, severall .authors:^'" have' investigated the effects of 
permeability variation. 

To the knowledge of authors, the effects of time dependent permeability variation 
along with time dependent suction velocity in unsteady flow have not been discussed, 
though these effects show their significant role in chemical engineering, 
turbomachinery and aerospace technology. Therefore, the aim of the present inves- 
tigation is to study the effects of permeability variation and oscillatory suction veloc- 
ity on free convection and mass transfer flow of a viscous fluid past an infinite vertical 
porous plate. The permeability of the porous medium is considered as K(t') = 
K'oCl+se'" ' ) and the suction velocity is assumed as v(t') = -Vo(l+ sa'” ') where Vq > 0 
and s « 1 is a positive constant. 


Mathematical Analysis 

We consider unsteady free convection and mass transfer flow of an incom- 
pressible, electrically conducting, viscous fluid past an infinite vertical porous plate in 
a porous medium. In cartesian coordinate system, :c'-axis is taken along the plate in the 
direction of the flow and y-axis normal to it. A uniform magnetic field is introduced 
normal to the direction of flow. In the analysis, we assume that the magnetic 
Reynold's number is much less then unity so that the induced magnetic field is 
neglected in comparison to the applied magnetic field. Further, all the fluid properties 
are assumed constant except that of the influence of density variation with 
temperature. Therefore the basic flow in the medium is entirely due to buoyancy force 
caused by temperature difference between the wall and the medium. Initially at t' < 0, 
the plate as well as fluid are assumed to be at same temperature and the concentration 
of species is very low. When t' > 0, the temperature of the plate is instantaneously 
raised to Tl and the concentration of species is raised to C^. By use of the 
Boussineseq's approximation, the governing equations in dimensionless form are as 
follows. 


Momentum equation : 


4 dt 



■¥ze 



G,r + G„C + 




M^u 


( 1 ) 
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Energy equation : 


4 dt 


(> 


+ se 



1 

P, dy^ 


Concentration equation : 


4 dt 




dc 1 ax 


dy 


Sc dy'^ 


( 2 ) 


( 3 ) 


where the symbols have their usual meaning as defined a Achaiya et a/.” 

The non-dimensional quantities introduced in the above equations are defined as : 


,_Voy .VoV 45n' «• __ c-c„ 

9 4a Vo' Vo’ c,-c„’ 


(^w — (Grashof number) Sc = —, (Schmidt number), 


Vn D 


is:. 


-,(parandtl number), M = — 


Vn 


(Magnetic parameter) 


and G„ = — — -'*'j — (Modified Grashof number) 

Vo 

The relevant boundary conditions in dimensionless form are : 

M = 0, r = 1+ee'"', C = 1+ sa'"‘, at y = 0 

M->0, r -> 0, C-^O, at y = 00 

In order to solve the system of equations (1), (2) and (3) under the boundary 
conditions (4) we assume 
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wCv,0=“oCP)+ewi(p)e'"' 

T(y,t)=niy)+zT,{y)e‘’'‘ 

and C(y,t)=Co{y)+zC^{y)e‘'" (5) 

Substitution (5) into the equations (1) to (3) and equating the harmonic and non- 
harmonic terms, we get : 

Uq+u'q- a,Wo = -G^Tq - G„Cq (6) 

uf +u[ -i32«i = -G^Ti -G„C, -u'q jp~“o (7) 

-^0 

T*+PX = 0 (8) 


T^^+PX-jPrn=-Prn 

Ct+SX'o=0 


(9) 

( 10 ) 


Cf+C,C{-^5,C,=-5,C; (11) 

4 

where the primes denote differentiation with respect to>'. 

The boundary conditions (4) reduce to : 

Mo = Ml = 0, To = Ti = l, Co = Ci = l, aty = 0 

mq ~ Ml — > 0, To — T\ — > 0, Co ~ C\ — > 0, as y — oo (12) 

Substituting the solutions of equations (6)-(ll) under the boundary conditions 
(12), we obtain ; 
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T(y,t) = 6-”'^+ e 


^1 

1 

e-«ny 

+ i^Pr.e-Pry 

LV « ) 


n 


(13) 


C{y,t) = e-^^^+e 


K n J n 



(14) 


w(3'> t) = (^3 + ^4 

+ e[a,se-'"^^ +a,3e-'’'-^ +0,46-^-^]^"’' (15) 

With convention that the real parts of complex quantities have physical 
significance in the problem, the velocity, temperature and concentration field can be 
expressed in fluctuating parts as : 

wCVjO ~ MoCy) + s(Mr cos nt-M, sin nt) (16) 

T(y,t) = 7’o(y) + e(Kr cos nt-K, sin nt) (17) 

C(y,t) = Co(y) + e(Lr cos nt-L, sin nt) (18) 

where 


Mr = (Nil cos ^37 + Ni2 sin B2y)e'’^^ -(Ni cos Biy + V2 sin 5 ,y)e'^i^ 
+ (Ws cos B^y + sin B^y) 6'^^+ N5 e"”^ + W7 V? 


M. = (Ni 2 cos Biy + V„ sin B^y) -(N 2 cos Biy + Ni sin 5,y)e''^>’’ 


+ (N4 cos B2y + Ni sin Biy) + Ne e""^ + Ns +N10 

Kr= { cos 5 ,y - -^^sin B^y 

\ n J 
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nt 


K, = ^sin B,>’ --^^cos 


Lr = I cos B^y - ^^sin 
n 


,-^y 


( AS ^ 

L, = sin Bjy ^cos B^y 

V « 




Hence expression for transient velocity, temperature and concentration field for 
% 

= —are : 


« y. 


n 

2n 




( 19 ) 


T y, 


n 

2n) 


= Toiy)-eK, 


y,— \ = Co(y)-eL, 


( 20 ) 

( 21 ) 


Skin-Friction, Rate of Heat and Mass Transfer 
Skin-friction coefficient (i) at the plate, in terms of amplitude and phase, is 


T = 




= Iq+ e I I cos(«/ + a) 


( 22 ) 


y=0 


Heat transfer coefficient (Nu) at the plate terms of amplitude and phase 








= P,+ e I I cos(nt + P) 


( 23 ) 
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Mass transfer coefficient (S*) at the plate in terms of amplitude and phase 








>.=0 


= 5'^+€|0|cos(«t + y) 


(24) 


Discussion and Conclusions 

In order to get physical insight into the problem, velocity field, temperature and 
concentration field, skin-friction coefficient, rate of heat transfer and rate of mass 
transfer have been discussed by assigning numerical values to various parameters 
encountered into the equations. To be realistic, the values of Prandtl number (Pr) are 
chosen to be P,. = 1 1.4 which represents water at 4°C, Pr = 0.71 and Pr = 7.0 which 
represent air and water at 20°C respectively. The values of Schmidt number Sc are 
taken for hydrogen {Sc = 0.22), oxygen (Sc = 0.66) and ammonia (Sc = 0.78 ). In these 
numerical calculations Grashof number for heat transfer Gr = 5.0, 10.0, 20.0, modified 
Grashof number for mass transfer Gm = 5.0, 10.0, 15.0, magnetic parameter Af = 0.5, 
1.0, permeability parameter Kq = 10.0, 20.0 and frequency parameter « = 5.0, 10.0, 
15.0 are chosen arbitrarily. 

It is of immense important to mention here that the series solutions obtained for 
velocity u, temperature T and concentration C are expended in powers of s. 

Therefore to judge the convergence of these series solutions, a computer program 
was setup to calculate +Kf , and + 1}, , under the conditions 

+M^ « uq yjKj +Kf « To and yjL] + 1}, «Co. In this manner it was 
observed that the solutions in series form are valid the above mentioned restriction. 

Fig. 1 shows transient velocity field at s = 0.002, Pr = 0.71 and nt = Jt/2. It 
represents the curves due to change in the numerical values of Grashof number for 
heat transfer Gr, modified Grashof number for mass transfer Gm, Schmidt number Sc, 
magnetic parameter M and permeability parameter Kq. Fig. 2 shows transient 
temperature and concentration fields at s = 0.002 and nt - Ji/2. It represents the 
curves due to change in Prandtl number Pr and Schmidt number Sc respectively. The 
numerical values skin-friction coefficient (i) in terms of amplitude |A/1 and phase tan a 
are tabulated in Table 1 for variations in Gr, Gm, Sc, M, Kq and n respectively. Table 2 
represents the numerical values due to variations in the numerical values of Pr and 
frequency n on the rate of heat transfer in terms of amplitude |i?| and phase tan p. 
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Table 3 represents the numerical values due to variations in the numerical values of Sc 
and frequency n on the rate of mass transfer in terms of amplitude |2| and phase tan y. 

Table 1- Values of amplitude \N\ and phase tan a for skin-friction coefficient 


Gr 

Om 

S’. 

M 

Kc 

n 

|A^ 

tan a 

5.00 

5.00 

0 22 

05 

10.0 

5.00 

30.00448 

-0.2545! 

10 0 

5 00 

0 22 

05 

10.0 

5 00 

49 46438 

-0.29118 

5.00 

10.0 

0.22 

05 

10.0 

5 00 ' 

41.16401 

-0.21461 

5.00 

5 00 

0 60 

0.5 

10.0 

5 00 

24.98689 

-0.22516 

5 00 

5 00 

0 22 

1 0 

10.0 

5.00 

01 66431 

1 30813 

5 00 

5 00 

0 22 

05 

20.0 

5 00 

47 36441 

-0 27686 

5 00 

5 00 

0 22 

05 

100 

10.0 

25 94993 

-0.40371 


Table 2 - Values of amplitude \R\ and phase tan p for rate of heat transfer. 


S. No. 

Pr 

n 

m 

tan p 

1 

0 71 

5.00 

01 48214 

0 29255 

2 

7.00 

5 00 

13 82601 

0.00671 

3 

114 

5 00 

22.67623 

0 00179 

4 

0.71 

10.0 

01.75014 

046111 

5 

0.71 

15 0 

02.00082 

0 55456 


Table 3- Values of amplitude |g| and phase tan y for rate of mass transfer 


S. No. 

S. 

n 

lei 

tany 

1 

0.22 

5.00 

0.64138 

0 58794 

2 

0.66 

5.00 

1.42431 

0.37536 

3 

0.78 

5.00 

2 04135 

0 31725 

4 

0 22 

10 0 

0.85512 

0 74467 

5 

0.22 

15.0 

1.02098 

0.80761 
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Fig. 2- Transient temperature and concentration fields. 


The conclusions of the study are : 
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(i) An increase in G, or G„ or ATo decreases the transient velocity field. 

(ii) An increase in Sc or M increases the transient velocity field. 

(iii) Comparison of transient velocity distribution curves show that the curves 
rise gradually after attaining minimum value near the plate. 

(iv) An increase in Pr decreases the transient temperature field indicating that 
transient temperature field falls more rapidly for water in comparison to air. 

(v) An increase in Sc decreases the transient concentration field indicating that 
transient concentration field falls gradually for hydrogen in comparison to 
oxygen. 

(vi) An increase in Gr or or Ko leads to an increase in the value of amplitude 
|iV| where as an increase in Sc or Af or n leads to a decrease in the value of 
\N\. 

(vii) There is always a phase lag in the case of skin-motion except for magnetic 
parameter M where phase lead is observed. 

(viii) An increase in Pr or n leads to an increase in the value of amplitude |i?|. 

(ix) There is always a phase lead in the rate of heat transfer. 

(x) An increase in Sc or n leads to an increase in the value of amplitude \Q\ . 

(xi) There is always a phase lead in the rate of mass transfer. 
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m 2 


P.+^Jp^+inP, =^i+/Bi. 

= 1 Sc+^|Sc +inSc =A2 +iB2 , 2m3 =l + ^l+4ai , 


1 r , w2 , ^ 

Mi, ='i[l + ^l+4a2]=^3+'^3, 02=^1+ — , 






“3 = — ^ ^ ^4 = ~T' 

/=?-P^+a, ’ Sc-Sc+O]’ 


/:o ’ ■ 4 

^5 = G;. + /C| ^ 


^6 +'C2. ^7 =(«3 +^4)f'«3 ’ 


a I + /C 2 aio = '^l+'^2= 2 J . 


flu = ^3 +iiV4 


2 ~, "12 

m 2 -m 2 -02 

“8 


= Ni + lN 6 =-:^ 


ai 2 = ^ 7 +iN^ =—z — ai 4 = Ng + iNiq = - 

p}-Pr-a 2 ' Si-S^-ai 

4P G 

ai 5 = Nil +/-^12 =<^10 ~^12 *"^13 Cl = 

* w 


C2=-^:;^£^ C3=a3 

^ n ’ 


Pr 

0 


C4=04|^-S, 


= Q^2 +Qr^l ^ _ Q^l ~^r^2 ^ _ ^2^4 +*Jw^3 

‘ -^2 , 1/^2 I ^ 1^2 , tr2 


+ Ki ’ 


/:f + Ki ’ 


Ki+K^ ’ 


£ 2 ^^ 16.2^ 

Kf+K^ ’ 


l 6 Ki+n 


M 6 =- 


Anaj 


2 -^ 2 . ■'° 16J^2+„2> 


w 16C3^6-4 oCi 

N'j ”” '" 9 

16^1 


^9 


16C4jf7 -4«C2 


,, l 6 CiK 6 + 4 nC 2 
Ag 0 '"" 

l 6 Ki+n^ ’ 

JV ■_ 16C2/:7 +4«C4 

I 6 K 7 +W^ 


+ A ^3 - 7/5 - 7/7 ~ i ^ 9 , Ni 2 = 7/2 +-^4 "’-^6 ““-^8 ~^ 10 » 
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a:, =Af -Bf-A] -a,, 




Kj =A2-BI-A2 -aj. 


K 4 =2A2B2~B2-~ 


^5 = — ffj3 -aj. K(, — — Pf —a\, — Sg —S q —Q] 


A, =^ + 


2 2 V 2 " 


^2 = Y + ^^y ^[sf^ + S, 


A2 =- + 

■5 O 


1 1 


2 2V2 


^ -] 

■^(1 + 4 ^ 1 )^ +«^ +(i+4£ii)r 


^j.U^7^-pV\ b2=^MM^-Sc 


^3=^^ y(l + 4aiL +K^ -(l + 4ai 
To =<23^^ +a45<; -'”3(^3 +^4). 


V = V;. +/V, = 



R = + iR^ = • 



tanP = - 


e=er+*e, =• 


tanY = - 


Nr = A^N I - B]N2 + ^2^3 “ B2N4 

-(^3^11 -S3Vi2)-m3V5-E^V7 -S^V9 

V, = V 2 + B] Vj + /42-?^4 + ^2^3 

-{A 2 N 12 + B2Nii) — m2N(, - E;-Vg -S^Vio 


n 


n n 


Qr - Ai'^ 


4S,B2 


. ^ r. ^ScAo 4S; 

and Qf = B 2 ^ + 

n n 




PROC. NAT. ACAD. SCI. INDIA, 74 (A), II, 2004 


Effects of time dependent suction on heat and mass 
transfer flow of a viscous fluid in rotating system 

ATUL KUMAR SINGH* and N. P. SINGH 

department of Mathematics, V S. S. D. College, Kanpur-208002, India. 

Department of Mathematics, C. L. Jain (P.G) College, Firozabad-283203, India 
Received March 14, 2002, Revised July 6, 2002, Re-Revised December 9, 2002, Accepted July 15, 2003 

Abstract 

In the present paper free convection and mass transfer flow of an incompressible, viscous liquid 
through a porous medium in a rotating system past an infinite vertical porous plate is studied 
In the analysis, the heat flux is considered constant and the suction velocity is assumed to be 
time dependent. Two different cases namely externally cooled plate {Gi, > 0) and externally 
heated plate (Gr < 0) have been analysed The solutions for steady primary velocity, steady 
secondary velocity, time dependent primary velocity, time dependent secondary velocity, 
temperature distribution and concentration field have been obtained using perturbation 
technique The results obtained have been discussed 

(Keywords • time dependent suction/mass transfer/rotating system/porous medium/heat flux) 

Introduction 

In recent years, free convection and mass transfer flow of viscous fluids through 
porous medium have attracted the attention of a number of authors in view of its 
application to geophysics, astrophysics, meteorology, physiology, aerodynamics, 
boundary layer control and so on. In addition, convective flow through porous 
medium have applications in the field of chemical engineering for filtration and 
purification processes; in the field of agriculture engineering to study the underground 
water resources; in petroleum industry to study the movement of natural gases oil and 
water through the oil channels / reservoirs; in metal industry to study the flow pattern 
of the melted solids at different temperatures. The development of aeronautical, 
chemical, mechanical and computer engineering during last few decades, on the one 
hand and finding oil wells in oceans on the other hand, have added stimuli to the study 
of viscous fluid flows. 

Ahmadi and Manvi* derived the equations of motion for viscous flow through a 
rigid porous medium. Gersten and Gross^ applied these equations to study the effects 
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of transverse sinusoidal suction velocity on flow and heat transfer along an infinite 
porous wall and Ram and Mishra^ used them to study the unsteady flow of a viscous 
fluid through porous medium. Raptis"* investigated the unsteady flow through a porous 
medium bounded by an infinite porous plate subject to a variable plate temperature 
and constant suction velocity. Raptis and Perdikis^ further studied oscillatory flow of a 
viscous fluid through porous medium by the presence of free convective flow. Several 
authors®'*^ have studied free convection / mass transfer flow in rotating system under 
various boundary conditions. Recently, Singh et al}^ have analysed MHD free 
convection and mass transfer in dusty flow of a rotating viscous liquid. More recently, 
Singh et al have studied free convection effect in the unsteady MHD flow of an 
incompressible viscous fluid in a rotating system past an isothermal vertical plate 
under constant suction velocity. 

The present paper deals with free convection and mass transfer flow through a 
porous medium in a rotating system past an infinite vertical porous plate with time 
dependent suction and constant heat flux. The effects of various parameters for two 
different cases namely externally cooled plate ( Gr > 0 ) and externally heated plate 
( Gr < 0 ) have been discussed with the help of graphs and tables. 


Formulation of the Problem 

We consider the unsteady flow of an incompressible, viscous liquid by the 
presence of free convection and mass transfer through a porous medium past an 
infinite vertical porous plate. A cartesian coordinate system (x, y, z) rotating uniformly 
with the liquid in a rigid state of rotation is considered and a constant angular velocity 
(0,0, fl) is taken into account about z-axis. The vertical porous plate is considered in 
the plane and z-axis is chosen normal to the plate pointing towards the flow medium. 
Besides, in the analysis a time dependent section velocity w = — wo[l+se’”*] has been 
taken into account normal to the plate (wo and n are positive constants). 

Therefore, for the present configuration, the momentum equations, diffusion 
equation and energy equation governing the flow, in non-dimensional from (ignoring 
the stars introduced ahead), are : 


du 

dt 


2£;v+'(l 


-l-se' 


’dz 


= G,r + G„C + 




( 1 ) 
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— +2jEv + (l + se"”' 
dt ^ 



d\ 

dz^ 


1 

— V 

k 


( 2 ) 


dt 


(l + se-”' 



1 d^C 
Sc dz^ 


(3) 


and 




/’.(l + ee""' 


9z 


= 0 


(4) 


The non-dimensional boundary conditions relevant to the present problem are : 

t>0; fF=0, — = -L, C=l+eZ 2 e‘”' at z =0 

dz 

JV ~ 0, T — ► 0, C — ^ 0, as z — ^ oo, (5) 


where u and v are the velocities in x and y directions, T is the temperature of the fluid 
r« is the temperature of the fluid far away from the plate, C is the concentration of 
species, Co is the concentration of species far away from the plate, p* is the 
volumetric coefficient of thermal expansion, p is the volumetric coefficient of thermal 
expansion with concentration, K' is the thermal conductivity, Cp is the specific heat at 
constant pressure, D is the chemical molecular diffiisivity, k is the permeability of the 
medium, q is constant heat flux and the other symbols have their usual meaning. 

The non dimensional quantities introduced in the above equations are defined as : 


zwl 


!i = „. 
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where G= ^ , (Grashof Number), 5"^ (Schmidt Number), 

K'v* ^ 


Du i^C„ 

£ = —, (Rotation Parameter), = — ^,(Prandtl Number) 

Vo K 


G„ = (Modified Grashof Number), L■^ = — ~ — , 

Vq T'w - ^00 


and L,= -^ — , 

C-C 

W '^00 


Assuming W = u + iv, the equations (1) and (2) transform to : 


■ + 2iEv - (l + se"”' )— = G,T + G„C + ^ --W 
^ ' dz ' dz^ k 


Solution of the Problem 

To obtain the solution of (6), (3) and (4), we assume 

W= Wi(z) + sW 2 (z)e-"‘+ 

T ^ T,(z) + sr 2 (z)e-"'+ (7) 

C = Ciiz) + sC 2 (z)e'"'+ 

Using (7) is equations (3), (4) and (6) and comparing the coefficients of harmonic 
and non-harmonic terms, we obtain six equations. The solution of these equations 
under the transformed boundary conditions yield : 

C,(z) = e'^"‘' (8) 
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where 


C,i2) = 




-K, S, -Sz 
e ^ H — -e ^ 


V 


Z/1 z 
T,iz) = -^e '■ 


U2) = ^e'"^' 

A 0 


, -P 2 -K^2 \ ( -^.Z 

W^{z) = KAe -e e ^ -e 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


1 ^ 2 ( 2 ) = ^, 2 ^'^'" 


+ Kise''’^" +K,^e'^^" +K„e 




16 *= 


(13) 


k 2 2 


l + J\-4L. 1 + Jl - 4(Z3 - «) 

^3 =4 +/ 5 , =-^4 K,=A,+iB, =— ^ 


Ks=A,+iB.= , K.=A^+iBi=—^ — ^ 

^ P,{P}-P,-L,y Sl-S,-L, 


Kn=-G„ 


c 

L 

« 


K,=- 


_ LiG^P^ 


K, 


Kg= As+iBs=-K3(Ks+K6), 

Ku^A^+iBy=Sc(^-KA Kn=A3+iB3= 


Kxo^ A6+iB6=LiGr+K5Pr 


K, 


Kl-K,-^L,-n) 
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K\i— Ag'^iBg— 




Ku- A\Q+iB\(i- — — , ? 

Kl -K, -(I 3 -n) 


K\s- 






Sl-S,-iL,-n) 


Kn= Ais+iBi3=-(Kn+ Ku+ Ku+ K,s+ 


^[V(l + 4ao)'+64£' +(l + 4ao)]' 


5 .= ^[Vo + 4aoy +64£^ -(l + 4 ao)]^ «o =-p 


^2=- + ^^|^V(l-477 + 4ao)^ + 64jB^ +(l-4«4-4ao) ^ 


^ 2 = 7 ^ 7(l-4« + 4ao)^+64£^ -(l-4n + 4ao)]\ 




A - 

M^+4E^Pr M^+AE^Py ' Ml+AE^' 


M^+AE^ 


As=BxiBi+B,)-AMi+Ai), 
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As—L\Gr+ A^Pr 


Be=B3Pr 


At=S,-^-A^, Bj=-S^B^ 

In j 


Mt+4E^ 


M^+4E^' 


K,M, 

Ml+4E^ 


Ml+4E^ 


AiMs + B^M^ _ B^Ms - A^M^ 

Ml+Ml ’ MUmI 


_ A^Mj-2EB^ 
' M^+4E^ 


D _B^Mt+ 2EA(, 


Ay Mg -lEBy 

Ml+4E^ ’ 


_ B-jM^ + 2ElAi 

Ml + 


^ 13 - — (Ai+Ag+Aio+Aii+An), 


Bi3~(Bs+B9+Bio+Bi X+Bu), Mi =Pr{Pr ~ -^r + “o) > 


M 2 = 5 ,^- 5 ,+ao, Af 3 =i(:i^-^i+n + ao, M 4 = - /sCz + n + tto 


A /5 = aI -BI - A^ + n + a^, 


A /(5 — 2A\B\—B\—2E, 


A/, = Pj - + n + tto 


A/8=5;-S,+n + ao 


Substituting the values of C\(z), Ctiz), T\(z), T 2 (z), W\{z), and W 2 iz), in (7), we 
obtain 


C = +8 11,- Ig-*!" + g- 


Bl B\Pf — T —nt 

—e ' +8 -^e ^ -Le ' e 
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and 


r = ^5 ! e |+i::6| e ■ ^ -e 


+ s 


-K., -K. 


-S , 


K,^e +K,,e ^ +K,,e ^ ^ K,,e ^‘+K„e 


(16) 


Hence, primary velocity (u) is : 


u(z, t) = Ui+ SMje =AJe -e cos 5,z 1 + B^e sin J5,z 


“■i4 *^S 

+ A^e '■^cos5,z-e )-54e"'^'^ sin5jZ 


+ se" 


A^e +A^e +^io^ cos5,z 


-Ai 


Z Z 

Aj^e '' +^j 2 ^ ^ +Ai^e ^ C0SB2Z 


(17) 


The secondary velocity (v) is : 


^(T) t) = Vi+ evjg =B^ e - e cos ^jzl + sin B^z 


+ B^ e cos B^z-e ]-A^e '^^sin^jZ 


+ se~ 


B^e -^B^e +5joe sinfijZ 


""5 i4 

5,,e '■^+5,2^ ‘^+Ai^e cosSjZ 


( 18 ) 
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Skin-Friction, Heat and Mass Transfer 
The skin-friction coefficient of the fluid (x/) at the plate 

few] 

J.=o 


^Ks(K3-Pr)+K6{S-K,)-z{K^Kn+K2Ku+K,Ku 


+PrK,,+S,K^6+K,Kn)e-"' ( 19 ) 

Hence, skin-friction (Xp) due to primary velocity is : 
Xp=A3(A\-Pr)-BsB\+A4(Sc-Ai)+B4B\ 

—se ”\K\A 3 '^ K 2 Ag+ AjAio— BiBio 

+PrAll+ScAi2+A2Ai3-B2Bi3] ( 20 ) 

Hence, skin-friction (xj due to primary velocity is : 

'^i~B3(A \—P r)—A3B\+B4{Sc—A 1 )— ^4^1 
~ze AiBiq— B\A\g 

^PrBn+ScBn+B2Ax3+A2Bx3] ( 21 ) 

The rate of heat transfer i.e. the heat transfer coefficient in terms of Nusselt 
number (A(,) at the plate is 
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which satisfies the boundary condition. 

The rate of mass transfer i.e. the mass transfer coefficient in terms of Sherwood 
number Sh is : 


S, 


Sz J.=o 


( 23 ) 




( s 

sf 

II 

1 

1 

CO 

^1 

•^2 

L « . 

r 

/ n 


( 24 ) 


Table 1- Rate of mass transfer; (e 

= 0 002 , l2= 1.0 and « = 0 . 05 ) 




t 

= 0.22 

5 c = 0.30 

5 c = 0.60 

5 c =0 66 


5 c =0 78 

1 0 

- 0.22092 

-0 30117 

- 0 60229 

- 0 66252 


- 0 78298 

20 

- 0.22087 

- 0.30112 

-0 60218 

- 0 66240 


- 0.78283 

3.0 

- 0.22083 

- 0.30106 

- 0 60208 

- 0 66228 


- 0 78269 

4.0 

- 0.22079 

-0 30101 

- 0 60197 

-0 66217 


- 0.78256 

5.0 

- 0 22075 

- 0.30096 

- 0.60188 

- 0.66206 


- 0 78244 

Table 2 - Skin friction due to primary velocity; (Cooling case at £ = 0 002 ) 

k 

E n 

Pr 

Sc G, 

Om 

t 


10.0 

1.0 0 10 

0.71 

0.66 10.0 

8.00 

1 0 

- 33 2582 

20.0 

1.0 0.10 

0.71 

0.66 10.0 

8.00 

1.0 

- 32 1045 

10.0 

2.0 0 10 

0.71 

0.66 10.0 

8.00 

1 0 

- 32.4167 

10.0 

1 0 0.05 

0.71 

0,66 10 0 

8 00 

1.0 

- 33 3728 

10.0 

1.0 0 10 

7.00 

0.66 10.0 

8.00 

1 0 

- 12 2919 

10 0 

1.0 0.10 

0.71 

0.78 10.0 

8.00 

1 0 

- 33 0048 

1 0.0 

1.0 0.10 

0.71 

0.66 15.0 

8.00 

I 0 

-43 8016 

10.0 

1 0 0.10 

0 71 

0.66 10.0 

120 

1.0 

- 39.3439 

10.0 

1.0 0.10 

0.71 

0 66 10.0 

8 00 

2.0 

- 33 2494 
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Table 3- Skin friction due to secondary velocity, (Cooling case at s = 0 002) 


k 

E 

n 

Pr 

5. 

Gr 

Om 

t 

^5 

100 

10 

0.10 

0.71 

0 66 

10 0 

8 00 

1 0 

5 7471 

20.0 

1.0 

0.10 

0.71 

0.66 

10 0 

8 00 

1 0 

6 8939 

10 0 

2.0 

0.10 

0 71 

0.66 

100 

8 00 

10 

8 4749 

10 0 

10 

0 05 

0.71 

0.66 

100 

8 00 

1 0 

6 3176 

10 0 

1 0 

0 10 

7.00 

0 66 

100 

8.00 

1 0 

9 7743 

10 0 

1.0 

0.10 

071 

0 78 

100 

8 00 

1 0 

6 4867 

10 0 

10 

0 10 

071 

0 66 

15.0 

8 00 

1 0 

3 7729 

10.0 

1.0 

0.10 

0.71 

0 66 

100 

120 

1 0 

10.5947 

10.0 

1.0 

0.10 

071 

0,66 

10 0 

8 00 

20 

5 7387 


Table 4- Skin friction due to primary velocity; (Heating case at e = 0 002) 


k 

E 

n 

Pr 

Sr 

Or 


t 


10.0 

1.0 

0.10 

0.71 

0 66 

-10.0 

-8 00 

1 0 

33 2582 

20.0 

1.0 

0 10 

0 71 

0 66 

-10 0 

-8 00 

1 0 

32 1045 

10.0 

2.0 

0 10 

071 

066 

-100 

-8 00 

1 0 

32 4167 

lO.O 

LO 

0 05 

0 71 

0.66 

-lO.O 

-8 00 

1 0 

33 3728 

10.0 

1.0 

0.10 

7.00 

0.66 

-100 

-8 00 

10 

12.2919 

10.0 

1.0 

0 10 

0 71 

0 78 

-100 

-8.00 

1 0 

33.0048 

10.0 

1.0 

0.10 

0.71 

0 66 

-15.0 

-8.00 

1.0 

43.8016 

10.0 

1.0 

0 10 

071 

0.66 

-10 0 

-12.0 

1.0 

39.3439 

10.0 

LO 

0.10 

0.71 

0 66 

-10.0 

-8.00 

2.0 

33.2494 
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Table 5- Skin friction due to secondary velocity, {Heating case at e = 0 002) 


k 

E 

n 

Pr 

5. 

Gr 

G„ 

t 


too 

1 0 

0 10 

071 

0 66 

-100 

-8 00 

1 0 

-5 7471 

20 0 

1 0 

0 10 

0 71 

0 66 

-10 0 

-8 00 

1 0 

- 6 8939 

10 0 

20 

0 10 

071 

0 66 

-100 

-8 00 

1 0 

- 8 4749 

10 0 

1 0 

0 05 

0 71 

0 66 

-100 

-8.00 

1.0 

-6 3176 

100 

10 

0.10 

7 00 

0 66 

-100 

-8 00 

1 0 

- 9 7743 

100 

1 0 

0 10 

071 

0.78 

-100 

-8 00 

1.0 

- 6 4867 

100 

1 0 

0 10 

071 

0 66 

-15.0 

-8.00 

1 0 

- 3 7729 

100 

10 

0 10 

0 71 

0 66 

-100 

-12 0 

1 0 

- 10.5947 

100 

10 

0 10 

0 71 

0 66 

-10 0 

-8 00 

20 

- 5 7387 


Discussion and Conclusions 

In order to get physical insight into the problem, the effects of permeability 
parameter (^), rotation parameter (£), Schmidt number (Sc), Grashof number (Gr) and 
modified Grashof number (G„) on primary velocity field and secondary velocity field 
in case of externally cooled plate (Gp-O) and externally heated plate (Gr<0) are 
studied taking numerical values of various parameters. To study the effect of Schmidt 
number (Sc) on primary and secondary velocities, the values of Schmidt number (Sc) 
are chosen for Oxygen (Sc = 0.66) and Ammonia (Sc = 0.78). The value of Prandtl 
number (Fr) is chosen as Fr = 0.71 which corresponds to air at 20°C and one 
atmospheric pressure. For concentration of species and mass transfer the values of 
Schmidt number (Sc) are chosen for Hydrogen (Sc = 0.22). Helium (S'c=0.30), Water- 
vapour (Sc := 0.60), Oxygen (Sc = 0.66) and Ammonia (Sc := 0.78). The numerical 
values of the remaining parameters are chosen arbitrarily. 

Fig. 1 shows primary velocity profiles (u) for heating case (Gr < 0) and cooling 
case (Gr > 0) at Fr = 0.71 to observe the effects of permeability parameter (k), rotation 
parameter (E), Schmidt number (Sc), Grashof number (Gr) and modified Grashof 
number (G„). It is observed that in case of externally heated plate the primary velocity 
(u) increases rapidly near the plate and after attaining a maximum value, it decreases 
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while in case of externally cooled plate the primary velocity (w) decreases near the 
plate and after attaining a minimum value it increases. It is also observed that an 
increase in Gr or results in an increase in the primary velocity field while an 
increase in k, E or Sc leads to a decrease in the primary velocity field in case of 
externally heated plate. The effects of these parameters in the case of externally 
cooled plate are just reverse. It is interesting to note that numerically, the effects of the 
parameters on primary velocity are equal in value but opposite in sign. 

It is worth while to mention here that the effects of permeability parameter (k), 
rotation parameter (E), Schmidt number {S^ Grashof number (Gr) and modified 
Grashof number (Gm) on secondary velocity profiles (v) for cooling case (Gr> 0) and 
heating case (Gr < 0) at Pr = 0 71 were observed on the screen of the computer and it 
was noted that in case of externally cooled plate the secondary velocity (v) increases 
rapidly in the vicinity of the plate and after attaining a maximum value, it decreases 
but in case of externally heated plate the secondary velocity (v) decreases in the 
vicinity of the plate and after attaining a minimum value it increases. It was also 
observed that an increase in k, Gr or G^ results in an increase in the secondary velocity 
field while an increase in E or Sc leads to a decrease in the secondary velocity field in 
case of externally heated plate. The effects of these parameters in the case of 
externally cooled plate were noted just reverse. 

The equation (14) shows concentration of species. It is obvious that for various 
values of Schmidt number {Sc) the concentration field decreases as z increases for 
given values of s, t, n and Li. It is also obvious that an increase in Sc results in a 
decrease in the concentration of species for all z. Thus concentration field is maximum 
in case of Hydrogen {Sc - 0.22) and minimum for Ammonia {Sc = 0.78). The 
concentration for the gases Helium {Sc = 0.30 ), Water-vapour {Sc = 0.60), Oxygen {Sc 
= 0.66) lies in the said range. 

Table 1 indicates the rate of mass transfer in terms of Sherwood number (5*) with 
respect to time (t) for Hydrogen {Sc = 022), Helium {Sc = 0.30), Water-vapour {Sc = 
0.60), Oxygen {Sc = 0.66) and Ammonia {Sc = 0.78). It is observed that an increase in 
Schmidt number {Sc) decreases the rate of mass transfer {S^) while an increase in time 
{() increases the rate of mass transfer {S/,). Thus the rate of mass transfer for Ammonia 
is highest in comparison to other said gases. 

Table 2 represents the skin-fnction {tpj due to primary velocity for cooling case at 
s = 0.002 to observe the effects of k, E, n, Pn Sc, Gr, Gm and t. It is observed that an 
increase in k, E, n, P„ Sc or t leads to an increase the skin, friction (tp) while an 
increase in Gr or Gm decreases the skin-friction (tp). 
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Fig. 1- Pnmary velocity profiles 
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Table 3 represents the skin-friction (t^) due to secondary velocity for cooling case 
to observe the effects of k, E, n, Pr, Sc, Gr, and t. It is observed that an increase in 
n, Gr or t leads to a decrease in the skin, friction (Ts) while an increase in k, E, Pr, Sc or 
G„ increases the skin-friction(Xs). 

Table 4 represents the skin-friction (Xp) due to primary velocity for heating case to 
observe the effects of k, E, n, Pr, Sc, Gr, G^ and t. It is observed that an increase in k, 
E, n, Pr, Sc or t leads to a decrease in the skin-friction (tp) while an increase in Gr or 
Gm increases the skin-friction (tp). 

Table 5 represents the skin-friction (-Cs ) due to secondary velocity for heating case 
to observe the effects of k, E, n, Pr, Sc, Gr, Gm and t. It is observed that an increase in 
n, Gr or t leads to an increase in the skin-friction (ts)while an increase in k, E, Pr, Sc or 
Gm decreases the skin-friction(ts). 
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Abstract 

Extreme vanations m atmospheric temperature results in the damage of exposed human skin 
and subcutaneous region The present paper deals with the study of mathematical and 
numerical computation of disturbed thermo-regulation with damage causing excessive fluid or 
water migration. In particular, the coetaneous region is assumed to be exposed to excessive 
heating for a long period of time which leads to bum injury in dermal skin The physical 
parameters like rate of blood mass flow, the rate of metabolic heat generation and thermal 
conductivity are taken to be variable while rate of sweat generation is a function of tissue 
temperature The mathematical mode! which consists of heat and mass diffusion equations has 
been solved using variational finite element method. 

(Keywords : human skin/thermo regulation/damage function) 

Introduction 

The exposure of skin region of human body to heat source either in direct contact 
or kept at a distance for a long period of time results in bum injury. The degree of 
injuiy depends both on the temperature of the source and time of exposure along with 
the thermal nature of the tissue. Mortiz and Henriques' demonstrated graphically the 
time and temperature required to produce a specific level of bum injury. Guyton^ 
pointed out that the average critical temperature at which the pain starts is 45“C. Stole^ 
has mentioned that the range of thermal sensation extends above and below the 
temperature levels of comfort into heat and pain on one side and cold and pain on the 
other side. Each limit is reached when the tissue is destroyed. For practical purposes 
these limits may be stated as instantaneous tissue temperature of 27°C as upper limit 
where the process of systematic progression of thermal-coagulation can be 
demonstrated'*’^ while the lower limit is mentioned as 0 to 4°C. 

In the present paper, we have estimated the damage by investigating the process 
of heat and water migration in dermal region by mathematical modelling in the form 
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of partial differential equations. For computation of the temperature distribution the 
use of finite element method has been made. The severity of bum has also been 
studied using the damage function defined in the following section. 

The Mathematical Model 

Energy equation including the effects of convection due to blood circulation, 
conduction and metabolic heat generation'^ 

div [k grad T] + mct {TA-T)+S=(i (1) 

in which k, m, c* and S are the thermal conductivity, blood mass flow rate, specific 
heat and metabolic heat generation rate respectively. Ta is the temperature at the 
deepest subcutaneous tissue. There is evaporation of sweat from the skin surface so 
that the equation (1) subject to the following boundary condition : 

K—=h(T-Ta)+LE (2) 

dx 

where, h is heat transfer coefficient, Ta is environmental temperature, and L and E are 
the latent heat of water (also sweat) and rate of sweat evaporation respectively. The 
bum injury causes (a) metabolic change (b) circulatoery change and (c) body fluid 
alteration. Due to bum injury large amount of fluid is lost from the body by 
evaporation through the bum wound. 

Due to water transport through tissue, eqn. (1) along with boundary condition (2) is 
then coupled with the following mass diffusion equation governing water magration. 

div [ D grad C] + :R=0 (3) 

in which D is the mass diffiisivity of water in the tissue, C is the water 
concentration in tissue and R is the rate of sweat generation. The appropriate 
boundary condition at the skin surface is given by 



D^=hu{{C-C,) 
ax 


( 4 ) 
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where ho is defined as the mass transfer coefficient and Ca is the water concentration 
environment. 

To complete the thermal (bum) injury model, a damage function devised by 
Henriques and Moritz^ has to be included. This damage function has subsequently 
been used quite widely by several workers^’*’’"’®. With this injury function the 
cumulative damage incurred during a bum can be predicted. 

The total injury, Q, at any point in the tissue is obtained by knowing the local 
transient temperature history, i.e., 


Q(x) = A Jexp 
0 


-AE 

RT(x,t) 


dt 


( 5 ) 


The function Q(x) was quantified by Henrique to identify various injury 
thresholds. At f2 = 0.53 is the starting value of irreversible epidermal injury, whereas 
at Q= 10'* complete transepidermal necrosis occurrs. In between there exist several 
stages of the thermal damage. The mathematical study of thermal thermo-regulation in 
in vitro tissue was taken up by several authors"'*^’ Saxena and his co-workers'*’^®, 
obtained interesting results for normal atmospheric conditions and abnormal growth 
in skin layers. Diller and Hayes^' modeled simple bum process using finite element 
technique in blood perfused skin. 

In the present work we have studied the thermal effect of hot atmospheric 
_ conditions on the SST region upto the extent of bum in one dimensional steady state 
case. The region under study has been assumed to be divided into three layers with 
total thickness c. The first layer, epidermis, is of thickness c measured from the skin 
surface; dermis, the second layer has thickness b and the third layer, the sub dermal 
tissue, is of thickness measured from the outer skin surface. The assumptions 
regarding external and internal conditions are given in Table 1. The boundary 
conditions are : 

At the skin surface (outer boundary) (x=0), we have the condition given by eqn. 
(2) and eqn. (4) with 

r= To 


and C =C 


( 6 ) 
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At the inner boundary (x = c) 

T= Tb (body core temperature) 

and Ta = Tb (7) 


Method of Solution 

Solving eqn. (3) along with boundary condition (4) in view of the values assigned 
to C, D and R in tissue layers and at interfaces conditions (Table 1), we obtain : 

C, =[l/2p(^)-a)/ A+(r, + T2)x + 1/2 \i(b-a)/hD(T^+T2) + Ca 

-\i,[\l2{,bT\^aT2)lb-a)x^+ \l'i)T2-T^)l{b-a)x^yD:^ 

■^[\l2]xlD^{b-a){b^TxH{b-2a)T2}]x 

+ \l2\i{b-a){alD,+\lhD) (T, + T 2 ) 

+ 1/6 iisilD2{b-a){i3ab-a^-3b^)T]+i6ab-3b^-2a^)T2)} 

Accordingly the rate of sweat evaporation E is estimated as 

£= ^lxib-a)iT,+T2) ( 8 ) 

where Ti and T 2 are unknown temperatures to be determined. 

Now coupling eqn. (1) with eqn. (2) and comparing with Euler's Largrange's 
equation, we get the following equivalent integral form ; 

I = +mcb{TA-Tf-2ST\dx+\l2h{T-Taf+ LET 

JL A oX 


( 9 ) 



Interior of the region At interfaces 
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The variational integral (9) for the whole region can be given as the sum of all h 
defined for each layer. Thus, 

/=!/, ( 10 ) 

»=1 


Evaluating the integral (9) for each layer by substituting layerwise expressions and 
interface conditions given in Table 1 and boundary conditions (7) and ( 8 ) respectively, 
we obtain: 


h = K^/ 2a(r, -Tof+l/2h{T(>-Taf+l/2Liiib-a)ToiTi+T2) 


l2=Fi(T2-Tif+F2ibT,-aT2)^+F,{T2-T,)+F,(bTi-aT2) 


+Fs(T^T0ibTi-aT2)+F6 


h=^E,(h-T2f+E2(cT2-bT2f+E2iT2-T2)+E,(cT2-bn) 


+Es{T2-T2){cT2-bT2)+F6 


The values of F/s , £',5 for /= 1 to 6 are given in Appendix I. 

Minimizing / for To, T\ and as = Tb, a constant, we get the following system 
of equations : 

XT=Y ( 11 ) 

Here 



A 

Xj 

13 


To 


W, 

X = 

M, 

M 2 

M 3 

; T = 

Tx 

Y = 

W2 


A^i 

N 2 

A 3 1 


T2 


W 2 


where L',s , M's and W’s W/s (/= 1,2,3) are defined in Appendix II. 
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Results and Discussion 


For numerical computation, following 
used 

K\ = 0.5 X 10'^ cal/cm sec°C; 

7 n 3 Ci= 0.525 X 10'^ cal/cm^ - sec°C; 
h = 3 x\0'^ cal/cm^ - sec°C; 

Tt, = 37°C; 

£),= 0.97 X 10‘^cmVsec; 
hd= 59.5 X 10"^ cm/sec, ; 


values of the quantities®’ have been 

jfiTa = 1.0 X 10'^ cal/cm sec°C 
S 3 = 0.3 X 10'^ cal/cm^ sec^C, 

I = 579 cal/g, 

R = 0.33 X lO''' gm/cm^/sec., 
£>2= 0.34 X 10’'' cmVsec, 

C„= 4.51498 X 10"’'*g/cm^ 


a = 0.10 cm; i = 0.35 cm; c = 0.50 cm, 

A = 3.x lO^Vsec.; AE = 6.3 x lOVKmol; R= 8.31 x 10^ J/Kmol°C 

Using these numerical values in set of eqn. (11), we get the values of nodal 
temperature and sweat evaporation at different values of heat source temperature 
(Table 2). 


Table 2- 


U°C) 

U°C) 

T/CQ 

T:(rC) 

E2 

(g/cm^-min) 

60 

50.97 

46.70 

41 06 

5.929 X 10-4 

65 

54.30 

49 07 

41.91 

6 147 X 10-4 

70 

57.63 

5143 

42.77 

6 365 X 10-4 

75 

60.95 

53.80 

43.63 

6.582 X 10-4 

80 

64.28 

56.16 

44.49 

6.8 X 10-4 


\ 

Nodal temperatures (T, i - 0.1.2.) and rate of sweat evaporation (E) at different values of source 
temperature (Ta) 
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Fig l~The curves numbered I, II, III, IV and V are drawn for Ta- 60°C, 65°C, 75°C and 80°C 
respectively Subscript 's' stands for skin surface and 'i' 
for interface of epidermis and dermis 

By substituting the nodal values of tissues temperature (Table 2) in 7^'^ (/=1,2,3) 
we get temperature profiles in dermal region at different values of source temperature. 
After the completion of steady state, we obtained damage function at skin surface and 
at interface of epidermis and dermis where nodal temperatures are above critical 
temperature. Using eqn. (6), Table 2 and corresponding numerical values of constants, 
we get damage ions(Fig. 1). 

From the Table 2 we observe that the tissue temperature T increases as the value 
of Ta increases and the rate of sweat evaporation also increases as value of Ta 
increases. We also find that the extent of bum increases as source temperature (Ta) 
increases. With increment in the extent of bums, the rate of sweat evaporation also 
increases. 

From Fig. 1, we observe that 

(i) At Ta = 60°C and at skin surface, irreversible burn is obtained after the time 
1 1 min. 24 sec and injury into entire epidermis may be available after the 
time 21 min. 24 sec. approximately. 
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(ii) At Ta = 60°C and at first interface (interface of epidermis and dermis), 
irreversible injury is available after the time 04 hrs 21 min. approximately. 
Injury into whole dermis may be available after the time 08 hrs 06 min. 
approximately. 

(iii) At Ta - 65°C and at skin surface, irreversible burn is obtained after the time 
1.04 min. and entire epidermis may be affected after the time 1.99 min. 
approximately. 

(iv) At Ta - 65°C and at first interface, irreversible injury is available after the 
time 45 min. and injury into whole dermis may be available after the time 
85.50 min. approximately. 

(v) At Ta = 70°C and at skin surface, irreversible injury is available after the 
time 6.20 sec. and entire epidermis may be affected after the time 11.6 sec. 
approximately. 

(vi) At Ta = 70°C and at the first interface, irreversible injury is obtained after 
the time 8.2 min. and entire dermis may be affected after the time 15.4 min. 
approximately. 

(vii) At Ta = 75°C entire epidermis may be affected within 1 .5 sec. and at the first 
interface irreversible burn is available after the time 1.5 min. and injury into 
entire dermis may be available after the time 2.85 min. approximately. 

(viii) At Ta - 80°C entire epidermis may be affected within 0.15 sec. and at the 
first interface irreversible burn is obtained after the time 1 7.2 sec. and injury 
into entire dermis may be available after the time 32.2 sec. approximately. 

From the Fig. 1, we also notice that initially the damages are linear, however, 

afterwards they show little fluctuations with time. 
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Appendix I 

Z, = [K\/a + h\\ Z2 = [112 p\iL(b-a)-KJla),. L3 = [172 \iL(b-a)] 
fVi = [bTaJ; W2 [Fs-bF,] 

W 3 =[( 2 E,+ 2 bcF 2 -fb+cJEs} Ti-F3+aF4+E3- c£4)] 

Ml = [1/2 ^lL fb - a)- Ki/aJ M2 = [Ki/a+ 2Fi + 2b^F2 - 2bF3],. 

M3 = [-2Fi - 2abF2 + (a^c) Fj] 

Ni = [1/2 jiZ [b-a)] 

N2 = [ -2F, - 2abF2 [a + b) F5] 

N3 = [2Fi + 2a^F2 -2aFs + 2Ei + 2c^E2 -2cE5l 

Appendix II 

Fi = [\IA(Ki +K3/(b-a) + \l2A(b-a) m3Cb(a^+2ab+ 3b^)\, 

F 2 = [\IA(b-a)m3Cb\-, 

F3 =[-\l\ 2(c-a) msCbTb (b-a) (a+3b) -1/6 S3 (a+2b)] ; 

F4 =[ -\l3(c-a) msCbTb (b-a) -1/2 Fj =[1/6 {b-a) msCb {a+2b)] 
Fe [ 1/8 (c-a)^ msCb T,2 (b-a)3\\ 

El =[l/2(c-b)K3 + l/6(c-b)m3Cb (b^+C^ +bc);. E2 =[ \l2(b-c) /WiC*],. 
E3 = [-l/6(c-a) m3CbTb {2b^ +2C^ + 2bc - 3ab - 3ac) -1/2 S3 (b-c)]; 
F 4 = [-m(c-a) m 3 CbTb {b+c-2a) - ^j]; Fj = [\l 2 (c-b)m 3 Cb{b+c)]-,. 

Ffi = [l/6(c-a) m3CbTb{b-a){a^ +b^ + 3c^+ab -3bc -3ac)] 
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Abstract 

A simple analytic proof of Mond-Shisha inequality for operators with a complete orthonormal 
set of eigenvectors is given. We then generalize the inequality for normal operators It is then 
generalized to a pair of selfadjoint commutative operators under certain conditions. The 
structure of the vectors for which one of the inequalities becomes an equality is also studied. 

(Keywords : strictly accretive operators/mean of order r/stationary vectors) 

Introduction 


Let q\, q2,—.qn be positive numbers with =1. For every sequence (xi, 

X 2 ,..Jc„), with all Xk> 0 and for every r 0, consider the mean of order r, Mr (xi, 

X 2 ,..jc„) defined as Jr* | • Let r, s be given reals r<s,rsj^0, then by a classical 

inequality (see Hardy et al.^ or Beckenbach and Bellman^) we have 


1 ^ Ms (Xi, X2,...Xn)/Mr (jT], X2,...X„). (1) 

Cargo and Shisha^ gave an upper bound for the ratio in (1). They assumed that 0 < 
a < p and a ^ ^ p, A: = 1, 2....n. Then setting y = p/a, they found 


Ml (JCi, X2,...x„)/Mr X2,...x„) < A 


( 2 ) 
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where A = 


f ’■if-f) 

n 


r 

1 

X, 

t 

1 

I 

[(r-. 



-l/r 


Mond and Shisha'* showed that these inequalities between means give rise to 
inequalities involving a bounded, selfadjoint operator^ on a Hilbert Space H. In fact, 
they proved the following. 

Theorem (Mond-Shisha) Let ^ be a bounded linear selfadjoint operator on a 
Hilbert space H with ml < A < MI mdO < m < M, I being the identity operator on H. 
Then 




( 3 ) 


where /e //, |1/[| = 1 and A is given by (2) which y = MIm. 

For other related results the reader may refer to Mond and Pecaric^"' ’ . 

In this paper we study the structure of the vectors that make {Ay f)^''‘l{A’'f,jf''' = 
(j)^(/) (say), for |j/[| = 1, stationary, which we call stationary vectors. The concept of 
stationary value and stationary vector was introduced by the first author Das’^ in 
studying the Rayleigh quotient of an operator and it was later used by us Das et al}^ in 
studying the anteigenvectors of a strictly accretive operator. We here give the 
definition of stationary value and stationary vector as follows. A function ^A{f) of 3 
unit vector /is said to have a stationary value at /if the function w/?) of a real variable 
t, defined by 


= 


f^tg\ 

_(a'(/+ tg),/+ tgf' [f + tg,f + tg)^'' 




has a stationary value at f — 0 for any arbitrary but fixed vector ge H. In other words 
M'g(0)=0 for all g. Then/is called a stationary vector. 

We show that a suitable linear combination of two eigenvectors is a stationary 
vector. 
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If ^ is a bounded linear selfadjoint operator with a complete orthonormal set of 
eigenvectors then we find the stationary values of A in terms of the eigenvalues of 
In case the maximum of is attained, the stationary vector is a linear combination 
of eigenvectors corresponding to the minimum and maximum eigenvalues. If the 
maximum is not attained, by our method we obtain an upper bound of ^A{f), which 
serves as an alternative proof of Mond-Shisha inequality for operators with a complete 
orthonormal set of eigenvectors. 

Next we generalize this inequality first to a strictly accretive {Re ^ > 0) normal 
operator and then to a pair of bounded, selfadoint and commutative operators. 

Throughout this paper H denotes a complex Hilbert space, with inner product (,) 
and norm || |1, 1 is the identity operator on H and .r4 is a bounded linear operator on H, 
and and e, denotes the eigenvalues and eigenvectors of A respectively. 

Stationary vectors of {A%ff'^l{Ay,f)^''' 

Here we find the necessary and sufficient condition for a unit vector / to be a 
stationary vector of ^Aif) when A is selfadjoint and show that for any two eigenvectors 
e\, ei there exists a suitable linear combination / = aiei + of e\ and ei such that / 
is a stationary vector of ^Aif)- Then we find the stationary values of (|)^(/) in terms of 
eigenvalues of A assuming that^I has a complete orthonormal set of eigenvectors. 

Theorem 1 : Let ^ be a bounded, linear, selfadjoint operator on H and e\, 62 be 
two eigenvectors of A with corresponding eigenvalues A.i, Xi. Then there exists a 
linearcombination / of e\ and ez such that /is a stationary vector of ^Aif)- 

Proof-. If (j)^(/) has a stationary value at a unit vector / then we have vt'g(0)= 0, 
where M'g(t) is defined as 


f^tg] 

(^" (/ + tg\ f + tgY {f + tgj-¥tg)\ ^ 


(A'if + tg),f 


For 11/11 = 1 set {Ay,f) = a" and {^ff) = b'', then w;(0)= 0 implies 


a^~'b 


Re(^V.g)- 


ab 


l-r 


>Re 



— a6Re(/,g) = 0. 
sr 


s 
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Since g is arbitrary we have 

rb' f - sa' A' f + (s-r) a'b''f= 0. 

This equation characterises the vector for which ^^if) is stationary. 
Suppose /= aiC] +a 2 e 2 , laip+ jaaf =1 
Then 


A^f= + ajA-jCj 

A''f= a,A,ie, +a2^2^2 
(Ay,f) = \afx\ +\a^\\\ 
(Ayf)= \a,\\\+\a,\X. 

Choose 


so that 


|a,|“- 


jX', - y, ) 




Now L.H.S. of (4) 


(4) 


= rb'Ay-sa‘Ay+ (s-r) if by 

= WlaipA.'j+iajpX^,) aai-s(|a,pAf,+|a2pXf2) 
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+ (s-r) (|aip Af,+|a2pA,2) (|ai|^X,'j+|a2pX2) cci} + 
{/•(laip A,'i+|a2|^A,2) a2A,2- s(|aipA,',+|a2|^A,2) a2A,2 
+ (s-r) (|aip Af,+|a2p Aj) (|aip)i^+|a2p Xj) a2} 62 

= 0. 

so that /is a stationary vector of ^A(j)- This completes the proof 
Also in this case 

(^y,/)= ia,pAl+la2fA.^2 

_ r(r2A,",-A.';A"2) 

t-r)^ -iQ' 

(A';,f)= |a,fX', +|ajfx, 

(s -r)^ -^2) 

so that 




f rl 

7 “') 

u 

n 

[ si 

k;"/! 

[(- 


T)J 





Next suppose A has a complete orthonormal set of eigenvectors ei with 
eigenvalues We now prove the following theorem. 


Theorem 2 : Let .<4 be a bounded, linear, selfadjoint operator on iif with a complete 
orthonormal set of eigenvectors ei and corresponding eigenvalues %i. Let r, 5, Mbe 
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given reals where d <m< M,r <s mdrs Let Ahe such that ml <A<MI. Then 
{Ay,f)'''XAj,f)^'' < A for/e //, 11/11 = 1, where A is given by eqn. (2) (y = MIm). 

Proof : Suppose / e i/, H/fl = 1 be a stationary vector of <t)^(/) and let/= 2(^6, )e„ Af 
= Substituting the values of f, Af'm (4) and equating to zero the coefficients 

of e, we obtain 

{rb^ X]-sa^ X[Hs-r)a^y) (f,e,) = 0 Vz 
If (f,e,) ^ 0 then 

rf X-'l+(s-r) = 0 (5) 

If /is an eigenvector then the characteristic eqn. (4) is satisfied for which a = b= 
X, where X is the corresponding eigenvalue. Then ^A(f) = 1> which gives the minimum 
value. 

If/ is not an eigenvector then the following cases arise. 

Case 1 : Eqn. (5) is satisfied by only two X,'s, in which case / is a linear 
combination of the corresponding two eigenvectors. 

Case 2 : Eqn. (5) is satisfied by more than two X,'s, in which case / is a linear 
combination more than two eigenvectors corresponding to two distinct eigenvalues. 

We shall however show that there always exists a suitable linear combination g, of 
two eigenvectors corresponding to two eigenvalues X, and Xj, X^Xj such that - 
(t)^(g) and is stationary at g. 

It suffices now to find a and a, such that 

g = a,e,+(Xjej, |a,p + |a^p = 1, 

= {Ay, f) = (A^g, g), y = (Ay f) = (A'-g, g) 

and (4) is satisfied for g. Now 

= ia,p r, + la/ and V = la/ AT, + |a/ . 
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Since |a,p+ |aj^= 1, we get 


, , b'-y:, 

1“'' " v/ 

A-i Kj Kj 

But since X, and Xj satisfy (5) we ha^^e 

ryyf,-sa^ (s-r) = 0, rb'X^j-sa^ X+{s-r)a’‘b'' =0. 


From these equations, it follows that 


a 


^ -X>') 


and b'' 


(s-r)(x'-X’J 


With these values it is seen that the two values of la,p are the same. Hence a 
suitable g is obtained. 

Thus we get ^A(f) in terms of the eigenvalues of .4. Next let 


/’(«) = 


f 

k. 

\ 

-"j 

1 ■ 

n 

f 



Ll 



K-i). 

L(^- 

-r)(«^-l)J 


where u = 


hL 


,X,>Xj. 


A little calculation shows thatF’fw) is an increasing function of u. 

So A (eqn. (2), where y = M/m) is an upper bound of ^A(f)- hi case the maximum 
of ^A(f) is attained, it is given by A and /= a,e, + afy, where e, and Cj are eigenvectors 
corresponding to the eigenvalues m and M resp. 
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This completes the proof. 

Extension of Mond-Shisha Inequality to Strictly Accretive Normal Operators 
Here we study the structure of the stationary vectors of 



where ^ is a strictly accretive (Re ^4 > 0) normal operator with a complete 
orthonormal set of eigenvectors. We obtain its stationary values in terms of 
eignevalues of A and taking the supremum of these values we get an upper bound of 

Theorem 3 : Let Abe a strictly accretive (Re ^ > 0) normal operator on H such 
that for any unit vector f s H, Af= 'L'k,{f,e,) e,. Let r and s be given reals where r<s, 
rs 0. If f is a stationary vector of (^a(J) then it is an eigenvector of ^4 or a linear 
combination of two eigenvectors corresponding to two conjugate eigenvalues or there 
exists a linear combination of two eigenvectors corresponding to two eigenvalues with 
unequal real parts and moduli, say g, such that 0 . 4 (g) = ^A(f) and O .4 is stationary at g. 
Further, a stationary value of ^Aif) is given by 

[ p.'’ cos'’ (j) - cos'' 0 ) ] [ 

[(s - r) (a,'” cos'' 0 - p'" cos'" <i))j I 

where and [xe"^ are the distinct eigenvalues referred to above and are such that 'k 
cos9 ^ p. cos(j). 

Proof : As in Theorem 1, 0.4(/) has a stationary value at a unit vector /if 
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Wg (t) = 


(j+tg),{/+tg) if+tgj+tgf"' 


^ A + A '^ 


l/r 


if + ^s\{f -^tg) 


has a stationary value at t - 0. For |j/f| = 1, set cf = ^■<Ja'A^ /,/ 


and f 


^A + A‘^ 
V 2 y 


/,/ 


Then we get 


ry (^^^ZAJ/■ 


sa 




/ + (s-r)aV/ = 0. 


( 6 ) 


This equation characterises the vector for which ^A{f) is stationary. Substituting 
Af= 'L'k,(f,e,)ei and A*f='E X , (f,ei)e, in (6) we obtain as before 


{ry\X,\^ - sy Re X]+ (s-r) yf} (f,e,) = 0 V*. 


If (fAi) ^ 0 then 

{ry |?i,f - sa" Re (s-r) yy) =0. (7) 

If/is an eigenvector then the characteristic eqn. (6) is satisfied for which a = |^|, b 
= Re X, where X is the corresponding eigenvalue. 

If /is not an eigenvector then the following cases arise. 

Case 1 : Eqn. (7) may be satisfied by only one pair (|(A.,|, Re 7-,), in which case 
{f,ei) 0 for two e,'s corresponding to two distinct X,'s. So the only way this can 
happen is that Af= X (f,e^ e, + X, (f,ej)ej. Here a = |A-1, = Re X,. 
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Case 2 : Eqn. (7) is satisfied by more than one pair (|(X,|, Re X,), then/is a linear 
combination of more than two eigenvectors corresponding to distinct eigenvalues. 


Proceeding as in Theorem 1 we shall however show that there always exists a 
suitable linear combination g, of two eigenvectors e, and ej corresponding to two 
eigenvalues with unequal real parts and moduli such that 0^(/)= and is 

stationary at g. 

Let X, - Xe'^ and Xj = pc'*, Xi^\L,X cos9 * p cos(j). It suffices now to find a, and ttj 
such that g = a, e, + afy, |a,p + |a;p = 1, rf* = f,f = g,g 




and (6) is satisfied forg. Now 


+lajpp*, b' = |a,pA,'' cos'O + |a,pp'' cos'’(j) 


Since |a,f + la^f = 1, we get 


2 j , I ,2 fc*" -p''cos''(l> 

l“'l -TT-^and also |a,p = ^ . 

-p. g cos 0 - p cos (|) 


But since X,, and Xj satisfy (7) we have 

ryx^ -sa^ X''cos'’0 + (s-r) a‘‘b''= 0, 
r6'’p^ -sa^ p'’cos'^(t) + (s-r) 0. 

From these equations it follows that 


p'’ cos'' (|) - p^X*" cos'' 6 ) ^ _ s(x‘p'' cos'" (}> - p^X'' cos'" 9 

{s- r) (x'' cos'" 9 - p'' cos'" (j)) ’ ( 5 -r)(x" -p*) 
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With these values it is seen that the two values of ]a,p are the same. Hence ag is 
obtained. Further alb gives a stationary value of ^>x(/). 

r-j — jC 

If the operator is positive then ^ = A and — ^ — = A also the eigenvalues are 

then real so that (|) = 0, 0 = 0 and Mond-Shisha inequality follows clearly as a 
corollaiy. 

Extension of Mond-Shisha Inequality to two Selfadjoint Commutative Operators 

Here we obtain the stationary values and the structure of the stationary vectors of 
{Ay, jf', where A and B are selfadjoint, commutative operators. The 

supremum of these stationaiy values gives an upper bound of the same. By taking 
P{A) and Q(A) in place of A and B respectively, where F(A) and Q(A) are polynomials 
in^, we have generalized it in yet another direction. 

Theorem 4 : Let A and B be two bounded, linear, selfadjoint, commutative 
operators on H with 0 < mj ^A< MJ and 0 < ^B ^ M^I. Let "k, and p, be the 
eigenvalues of A and B respectively with corresponding eigenvectors e,. For any unit 
vector/ e H, \eAAf= i:\{f,eM„ Bf= and {Ay,f)^'%TTf,ff‘' = H/(/) (say), 

where r and s are reals with r<s,rs=^0. 

If/is a stationary vector of v(/(/) then either it is an eigenvector of A and B both or 
a linear combination of two eigenvectors corresponding to two distinct pairs of 
eigenvalues of A and B, else there exists a linear combination of two eigenvectors, 
corresponding to two distinct pairs of eigenvalues, say g, such that v)/(g) = \\/{f) and 
is stationary at g. Also a stationary value of v|/(/) is 

where {X„ p,) and (kj, p^) are the two distinct pairs of eigenvalues of A and B referred 
to above. 

Proof : It can be shown that if / is a stationaiy vector of v|/(/) it satisfies the 
following equation 
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rf A’f - sa^Bj + (s-r) 0. (8) 

where cf = (AJ,/), f = (57, y) lL/11 = 1 • Substituting ^7= Z X] if, e,) e., and 57=Z (f, 

e,) e, in (8) we get 


rf }f-sa' \i'', +(s-A-) =0. 


( 9 ) 


if 7e,) 5 * 0. If/is an eigenvector then eqn. (8) is satisfied, for which a = X,andb = g„ 
where (k„ (i,) is the corresponding pair of eigenvalues. 

If/is not an eigenvector then the following cases arise. 

Case 1 : Eqn. (9) is satisfied by only two pairs of (K ^l,), in which case /is a 
linear combination of the corresponding eigenvectors. 

Case 2 : Eqn. (9) is satisfied by more than two pairs of (X„ p,)= then /is a linear 
combination of more than two eigenvectors corresponding to distinct pairs of 
eigenvalues. Proceeding as in the previous theorem we can show, that there always 
exists a suitable linear combination g, of eigenvectors corresponding to the distinct 
pairs (X„ p,) and (Xj, p,) such that v)/(g) = ^f(/) and ij; is stationary at g. Also the 
stationary value of vj/(^ can be calculated similarly. 

Corollary 1 : If B = A, r = -\, s = 1, the Kantorovich inequality follows (see 
Greub and Rheinboldt*^, and Kantorovich'^). 

Corollary 2 : If 5 = .^4, Mond-Shisha inequality follows. 

Corollary 3 ; By taking P(A) and Q(A) in place of A and B respectively in 
where F(A) and Q(A) are polynomials in A of the same degree, with positive 
coefficients so that P(A) > 0, Q(A) > 0, we get the stationary values of (5(^)7 
y)'7(g(y4)7y)'''' which is given by 




-l/r 


Thus we have extended the Mond-Shisha inequality in two directions. In Theorem 
3 we considered a strictly accretive normal operator A in place of a selfadjoint 
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operator A whereas in Theorem 4 we considered two selfadjoint operators A and B in 
place of a selfadjoint operator ^4. Since positive operator comes as a particular case of 
strictly accretive normal operator so Mond-Shisha inequality follows from Theorem 3 
and the same inequality follows from Theorem 4 by taking^ = 5 in Theorem 4. 
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Abstract 

Earth current variations related to earthquakes have been reported since the late 19^*^ century 
One of the earliest reports pertain to 1872 Bhavnagar earthquake which was preceded by strong 
currents flowing through telegraph lines. The other case involves the 1897 Great Assam 
Earthquake when the telegraph operators at Shillong experienced strong electric shocks prior to 
the incidence of aftershocks in the Meghalaya plateau. In that plateau, hundreds of millivolts of 
electromechanical potentials observed on hills showed significant changes prior to earthquakes 
in nineteen eighties The Russian scientists predicted many earthquakes m Kamchatka area 
analysing variations of earth currents and they also detected pressure induced polarization 
while rupturing large granite samples in the laboratory A programme based on the seismo- 
electrical signals predicted 18 out of 23 earthquakes in Greece. In China, the variations of 
geopotentials were used to predict the M 7.3 1975 Haicheng and many other earthquakes 
However, the natural ground potentials are vitiated by rains, extraterrestrial and man made 
cultural noises unlike the georesistivity methods where the noises can be filtered out 

The georesistivity investigations conducted by the geoscientists on the San Andreas fault, 
California and Garm region, Tadzhik republic since the nineteen seventies near the focal points 
using powerful deep resistivity units probing to depths of kilometers have shown that the rock 
resistivities decrease to a large extent before earthquakes. This has been corroborated by the 
results of rock fi*acturing experiments in America They lead to enunciation of dilatancy- 
diffusion models. 

The Chinese scientists have been monitoring the variation of earth resistivities covering the 
seismically active northern parts of their country by a network of more than 100 geoelectric 
stations. They have also conducted m situ experiments to relate minor changes in strains to the 
pattern of changes in resistivities of rocks. Based on the observed georesistivity precursors to 
more than 30 strong events over 25 years, they have developed travel time curves for the front 
of resistivity anomalies to locate the would be epicenters of the impending earthquakes They 
have deduced the magnitudes of the earthquakes in advance from three different parameters 


This paper is dedicated to the mem'ory of people who died in 26*** January, 2001 Bhuj Earthquake 
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such as the duration of the precursor, length of total anomalous area and maximum amplitude 
of the anomaly at the epicenter It has been possible to determine the medium term, short term 
and imminent stages of earthquakes from the trends of georesistivity anomalies. The imminent 
time of the earthquake has been determined from the appearance of the resistivity pulses 
induced by crustal tides of weakened crust on the verge of breaking supplemented by sharp 
changes in the self potentials 

At sensitive sites selected on geological and geophysical considerations in Shillong and 
Nagpur, the monitored earth resistivities have been observed to change and return to pre- 
anomaly level immediately before earthquakes occurring at 40 to 400 Ion from the sites The 
maximum change in magnitude of resistivities has been shown to be proportional to the 
hypocentral distance 

It is possible to integrate the methods developed by the Chinese and Japanese scientists with 
the indigenous methods and predict earthquakes in India using simple instruments. 

(Keywords . earthquakes/geoelectrical precursors) 


Introduction 

Earthquakes are generally caused by sudden movement within the lithosphere to 
release the strain that accumulates gradually over a period. The high pressure 
conditions and strains developing in the focal regions before and during earthquakes 
manifest themselves in changes of physical parameters of rocks and can be measured 
in the epicentral regions as well as at sites on ground surface several kilometers away 
from the epicenter. These parameters include seismic velocities, magnetism, self- 
potentials, electrical resistivities, radon concentrations, land deformation, tilt and 
strain and geohydrological conditions'. The changes in these parameters have been 
monitored to predict earthquakes by many geoscientists. In this paper, the discussions 
will be restricted to geoelectrical ( earth resistivities and self potentials) indications as 
they have been found to be more effective in forecasting earthquakes in China^ 
Japan^, Greece'', erstwhile USSR USA’ and also in India®. These studies are 
important for our country as six disastrous earthquakes have occurred during the 
thirteen years between 1988 and 2001. The earthquake striking Bihar Nepal border on 
20* August, 1988 killed 1000 people. This was followed by the Uttarakashi 
earthquake on 20* October, 1991 with a death toll of more than 1000 people. The 30* 
September 1993 Latur fcarthquake damaged nearly 4000 houses killing 9000 and 
injuring 20,000 people. The death tolls were 43 and 1 10 in 22"^' May, 1997 Jabalpur 
and 23"* April 1999 Chamoli earthquakes respectively. The 26* January, 2001 Bhuj 
earthquake demolished lacs of houses killing 195000 and injuring 17,000 people. It is 
true that loss of life and property can be minimiaed by insisting on construction of 
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earthquake resistant buildings, but it is not possible to replace or strengthen all the 
weak buildings in crowded cities located in earthquake prone belts of this poor 
country in near future. Hence, there is an urgent need to predict earthquakes. 

The seismic zoning map of India proposed on the basis of past earthquakes, 
seismicity and geological features is a process of forecasting the risk factors and 
indirectly earthquakes. The same process can be improved and supplemented by 
information regarding the ground conditions based on direct measurement of 
parameters to identify the vulnerable areas where earthquakes are likely to occur so 
that proper arrangements can be made in advance to mitigate the sufferings of the 
people. The objective of this paper is to describe some significant observations and 
relevant experiments carried out in the laboratory and in the field and highlight the 
fact that earthquakes can be predicted through scientific methods. 

In the self potential method, the potential difference is measured between two 
fixed electrodes using sensitive voltmeters or potentiometers to detect the electrical 
signals related to earthquakes. Unfortunately, these signals are vitiated by artificial or 
cultural noises such as leakage from power lines, and natural noises such as rainfall, 
solar flares and ionospheric currents. In contrast, resistivity methods employs an array 
of four or more electrodes. While current is passed into ground through a pair of 
buried electodes, the potential difference induced by electric field is simultaneously 
measured by other pair of electrodes. By varying the magnitude of current and 
fi’equency, the extraneous noise can be effectively eliminated. Hence, the earth 
resistivity methods are considered superior to the SP method. It is always desirable to 
work at sites away from the towns and combine these two methods as far as possible. 

Earthquake related Geoelectric Potentials 

Old classic cases : Yamazaki^ has summarized many classical and modem reports 
pertaining to relation between natural ground potentials and earthquakes. According to 
him, the Indian (Bhavnagar) earthquake of 14* December, 1872 was the first which 
was reported to induce strong currents through telegraph line connecting Valentia with 
London just before the event. Intensive studies on earth potentials were carried out in 
Japan by Noto in the nineteen thirties and subsequently by Yoshimatsu in the nineteen 
forties and fifties. They observed significant changes in earth potentials associated 
with earthquakes. The studies were followed up by the Earthquake Research Institute, 
Japan in the nineteen sixties. Rikitake and others buried spiral lead tubes in the ground 
at depths of about Im keeping the distances between the N-S and E-W pairs of 
electrodes 64m and 72m respectively. They observed a sudden change in earth 
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potentials for an earthquake of magnitude 5 centered at Matushiro on 23'^ November, 
1965, While the N-S pair of electrodes recorded 22 millivolts, the E-W pair recorded 
30 millivolts of earth potentials. 


SP studies in Kamchatka area 

In Kamchatka area, Fedotov et. al.^ observed earth-potentials as large as 
50mV/km several hours prior to earthquake as well during the events between 
September 1967 and February, 1968. According to them, anomalous geoelectric 
potentials might be related to stress variations that caused piezo-electric effects in 
rocks at the focal regions. In the same area, Sobolev® took continuous observations of 
electrotelluric field as a part of prediction studies. Potential differences were 
monitored at six stations along N-S and E-W directions maintaining a separation of 
200m between the pairs of electrodes. They observed changes varying from 8 to 13 
millivolts in geoelectric potential about one week before a deep earthquake of 
magnitude 6 at a distance of about 250 km on 27.05. 1972. Sobolev could predict 8 out 
of 12 earthquakes combining earth current observations with seismic wave velocities 
and statistical analysis. Later, Sobolev et. al. carried out experiments in laboratories 
rupturing large granite samples and observed that the samples emit transient electric 
current before rupture. Varotsos and Alexopoulos'* explain that the transient electric 
current comes from the orientation of the dipoles due to pressure induced polarization 
under changing stress on rocks at the focal regions. 


SP monitoring in Haicheng and Tangshan (China) 

An abnormal earth potential change was observed at a site about 20 km from 
Haicheng ( China) before the M 7.3 earthquake on 4* February, 1975. The equipment 
comprised two pairs of electrodes one of copper and other of graphite buried at a 
depth of about Im with a separation of 50m for E-W and 60 m for N-S electrodes. A 
group of geophysicists experienced in electrical prospecting monitored the potential 
differences using a self balancing potentiometer. As shown in Fig. la, the potential 
differences gradually increased to 170 mV for the N-S pair of electrodes. It is 
remarkable that the differences decreased to —90 mV the day before the earthquake 
(Fig. lb), and then suddenly increased to zero before the shock. The sharp changes 
observed on 2"^ and February, 1975 were accompanied by foreshocks and the 
earthquake was successfully predicted. 
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Fig. 1(a) Geodectrical potentials observed at dist. 20 km from Haicheng (China) from 5*^ January to 
5^^ February 1975 before M 7.3 4*^ February 1975 Haicheng earthquake; (b) The kinks on detailed 
graph on 2”^^ and Februaiy 1975 related to foreshocks. The sharp decrease of potentials 
started on 2”*^ February. This decrease and foreshocks helped in predicting 
earthquake (after Zhao et al , 1996). 
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Besides the work done by experienced geophysicists many amateurs participated 
in the programme. At every site, the amateurs buried two pairs of electrodes of 
dissimilar materials in N-S and E-W arrays as reported by the U.S Haicheng 
Earthquake Study Delegation. They monitored telluric currents using microammeters 
and helped the scientists. Fig 2 shows the self potentials observed at the station Xiji in 
the N-S and E-W channels The sharp changes and slow recoveries are attributed to 
the disposition of faults in the vicinity of the station in relation to the movement of 
fluid due to the M 7.8 July 1976 Tangshan earthquake^''. 


M 78 




Fig 2 (a) The self-potentials in N-S channel were observed at the station Xiji from January 1975 to 
April 1977. Mark the anomalous sharp change in S.P. before M 7.8 Tangshan (China) 
earthquake; (b) S.P . m the E-W channel at the same station showed 
decrease, (after Zhao et. ai, 1997). 
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Seismic electrical signals in Greece 

The changes in seismic electrical signal of the earth preceding earthquakes were 
studied in a network of 18 stations in Greece between 1981 and 1984 for earthquakes 
with magnitudes from 3.2 to 6.8'*. The electric field beyond the usual disturbances due 
to atmospheric and magnetic storms was observed measuring the potential differences 
between two pairs of electrodes in E-W and N-S directions at each station. A pair of 
brass ( or lead or non-polarized) electrodes were buried at depths of about 2m at 
separations varying from 30 to 200m. The stations were installed as far as possible 
from electric power sources to keep the cultural noise less than 0.2 mV. The potential 
differences observed at stations telemetrically connected to a center located at Glyfada 
about 15 km from Athens were measured after amplifying and filtering out 
frequencies higher than 0.3 Hz and depicted in multipen recorders or strip chart 
recorders with speeds varying 1-20 cm/h for each pair separately. 

The precursor seismic signals occurred 6-115 h before the earthquakes and had 
duration varying from 1-90 min. The signals appeared as transient changes of the 
potential differences between the electrode pairs measuring upto a few millivolts 
which are related to magnitude of earthquakes, epicentral distances and local 
inhomogeneities. Fig. 3 shows the seismic electric signal recorded at PIR-station at 
20-50 GMT on 17'*' February 1983. After about 56 h, an earthquake M4.8 occurred at 
a distance of about 10+5 km from the station. 

The current density in a particular direction was determined dividing the change 
of potential differences with a suitable effective resistivity. Measurements and 
calculations in two directions (E-W, N-S) gave the relative signal densities at that 
station. Considering that the signals are inversely proportional to the epicentral 
distances, the epicenter was located using the ’relative signal densities with an 
accuracy of 100 km. The magnitude of earthquake was deduced empirically with an 
accuracy of 0.6 units considering the relative signal densities and distances. Varotsos 
and Alexopoulos (1984)'* predicted 21 out of 23 events with magnitudes more than 5. 


Electromechanical potentials of Shillong Plateau 

Geoelectrical potentials were monitored over a hillock located in the reserve 
forest near Shillong from February to August, 1980. The site was selected on the basis 
of results of a systematic geophysical investigation carried out to detect occurrence of 
sulphide bearing rocks. High negative potentials were observed over quartzites 
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exposed on hill and the S.P profiles appeared like reflected images of the hills. The 
magnitude of potential in the vicinity of a peak point of a hill at a height of 1520m 
was - 1940 mV^^. The quartzites strike NE-SW and dip steeply towards south near the 
site. Since these rocks characterised by very high resistivities are devoid of graphite 
and sulphides as proved by drilling, the self potentials appeared to be caused by 
electromechanical properties of rocks". 
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The high electromechanical potentials were monitored keeping one potential 
electrode at the center of negative potential and the other electrode in a valley at a 
slope distance of 350m. As shown in tig. 4, the potential difference decreased from 
1500 mV on 2-3-1980 to 1340 mV on 9-3-1980. An earthquake of 4.5 occurred on 
9-3-80 at a distance of 55 km from the station, and after that the self-potentials 
increased. It is rather difficult to correlate the trend of anomaly with the instant of the 




MAGNITUDE ^ 


EARTHQUAKE 


epicentral DISX 

55 Km. 


Fig 4- Variation of apparent resistivity, self-potential, water level in a dugwell and rate of spring 
discharge before M 4 2 earthquake near Shillong (Nayak et al 1983) The earthquake 
occurred when resistivities approached pre-anomaly level and S.P. was low 
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earthquake, as measurement was taken only once daily. The potentials returned to the 
pre-anomaly level 4 days after that event. Another earthquake related potential change 
was observed for the earthquake that was felt in the area on 26-4-1980. No clear 
correlation was observed during the rainy season when the potential difference 
decreased to 900 mV. 

The electromechanical potentials are generally due to formation of a double layer 
of ions at the interface of the mineral particles and liquid. The potential difference is 
caused by displacement of the liquid phase of the double layer or movement of water 
through arenaceous rocks"’’^. It was envisaged that the height and steep slopes on the 
hill flanks provide hydraulic gradients for movement of water through permeable 
formations causing large potential difference. Although this phenomenon has been 
observed at other places, the magnitudes are much less than that at Shillong. This high 
earth potential may also be caused by electrification of rocks either due to piezo- 
electric effects’ associated with an earthquake or electro-elastic effects associated with 
an earthquake source'. It cannot be ruled out in this seismically active area where as 
many as 25-30 micro-earthquakes are recorded per day with about 2 medium events 
per month. In this context, it is worthwhile to quote from the classic work by 
Oldham'^ in his memoir on the Great Assam Earthquake of 12"’ June, 1897. 
“Telegraph workers experienced many electric shocks, some of them of considerable 
severity. Nearly all these shocks were preceded, but more often followed by an earth 
tremor, which rendered it impossible to utilize the earth as a return path for the 
current. The intensity of these currents was too great for them to be explained as a 
purely mechanical results of the disturbance of soil in contact with the earth plates. 
Though this is known to produce slight electric currents, the cause must be looked for 
at greater depth and more directly connected with the origin of the disturbances”. 

It is necessary to investigate the change in geoelectric potential related to the 
earthquake, whether it is entirely due to the charged rocks at the focus or due to 
changes of characteristics of rocks at the measuring station by taking simultaneous 
observations at a number of sites. If it is due to the former, the magnitude of change in 
potentials will be inversely proportional to the distance from the epicenter as in 
Greece. If it is due to the latter phenomenon, the magnitude of S.P. will depend on the 
rock type and hillocks irrespective of the distance from the hypocenter. 

Results of Experiments on Georesistivities under Controlled Conditions 

(1) Changes in electrical resistivities as a function of compressive stress in many 
crystalline rocks were observed by Brace and Orange'’*. The samples were subjected 



EARTHQUAKE PREDICTIONS USING GEOELECTRICAL PRECURSORS 


201 


to confining pressure upto 5 kb and to pore pressure of water of 500 bars. In majority 
of rocks, resistivity increased slightly upto about half the fracture stress. Beyond half 
and within 20 percent of the fracture stress, resistivity dropped sharply becoming one 
tenth of the initial value corresponding closely to an increase in porosity or dilatancy 
which took place under compressive stress. According to them, at low compressive 
stress, all cracks are closed resulting in decrease of the conducting paths, and hence 
the resistivity of rock samples increase. At stresses, prior to the fracture stresses, new 
cracks form in the compressed rocks and the electrical resistivity of rock sample 
decreases if conducting fluids fill these cracks. In another set of laboratory 
experiments Brace showed that the changes in rock resistivities in the focal region are 
entirely due to dilatant volume and not due to liquid to vapour transition of pore fluid 
as envisaged by some workers. 

(2) Scholz'^ showed through experimental study of fracturing process in brittle 
rocks that at about 92% of the fracture stress, a rapid acceleration of microfracturing 
occurs, and there is marked clustering of microfractures on the eventual fault plane. 
This suggests that there is increase of stress along the potential fault plane and 
relaxation at other points immediately before the onset of fracture and movement. It 
may help in indicating the time and place of imminent fracture. 

(3) Rikitake*^ observed a certain relationship between changes in earth resistivity 
and earth strain caused by tidal loading in 1952. Yamazaki^ undertook compression 
tests in the laboratory in 1965 which revealed that resistivity changes are larger than 
the mechanical strain by a factor several hundred particularly for porous and 
permeable rocks. He suspected that the contact section between the pores increases 
when the specimen is compressed giving rise to tremendous increase in the 
conductivity. 

(4) Zhao et. al}^ conducted in situ experiments in the field to determine the 
relationship between the strains of partially saturated soils to the variations of the 
resistivities. They applied horizontal pressure on walls of soil layers in trenches. Two 
hydraulic spiral jacks of 10-20 tons were used to exert pressure through steel plates on 
thick boards fixed to the walls of the trenches (Fig. 5a). To study resistivity changes of 
extended (bent) soil layers, two long parallel trenches were dug through soil layer and 
a bending moment was exerted on the soil body between the two trenches (Fig. 5b).To 
study the resistivity changes due to compression on rock layer, they selected a tunnel 
at depths of 100-200 m (Fig. 5c). The compressional stress was exerted vertically or 
horizontally by jacks on arenaceous limestones and electrodes were fixed on the side 
walls of the tunnels. 
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PRESSURE 


Fig. S-Schematic diagram of the experimental set-up using jack as a loader (the arrows represent the 
force source in the area (ANMB-the electric system), (a) Schematic diagram illustrating the 
loading of soil layer* and the vertical view (top) and sectional view (bottom) of 
the electrode arrangement, (b) Measurement points in the bending of soil 
layer experiment (c) The arrangement for measurement on the rock layer 
in a tunnel with cracks on side walls in the compressional experiment 
using jack as a loader (circle denotes the location of strain gauge) 

Fig (d) & (e) The resistivity on soil decreases with increase of 
pressure and vice-versa (after Zhao et al , 1990) 


To measure the strain, a rubber plug with an inner strain gauge was buried at 
depth of 22.5 cm at the center of the four electode system sticking the resistant strain 
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gauge on an intact part of the rock surface. The distance between the electrodes varied 
from 0.9 to 1.5m so that the deformation was uniform inside the volume of soil/rock 
to be tested. At each measuring point 3-4 profiles parallel, perpendicular and at 45° 
angles to the pressure direction (Fig 5a) were taken. 

Zhao et al}’’ observed that resistivity decreases (or increases) when the soil layer 
is compressed (or dilated) The relative resistivity variation is sensitive to strain 
changes, and smaller the strain, the higher is the amplification factor The resistivity 
change is directly proportional to intensity of force, size of the source, but inversely 
proportional to the distance between the source and measuring point The resistivity 
change shows a delay as this distance increases. It was also observed that the isotropic 
soil layers develop resistivity anisotropy due to stress. The resistivity decreases on 
contraction and increases with extension of the volume under testing (Fig. 5d & e) 

The experiments are very useful to understand the apparently anomalous 
phenomenon such as simultaneous increase and decrease of resistivities before 
earthquakes, appearance of anisotropy at the measuring point and propagation of 
anomalies from the epicenter. The experimental results have shown that the 
resistivities are very sensitive to changes in strain, and the area of the anomalous zone 
increases with increase of force at the focal region. It also explains the phenomenon 
observed in Japan by the Yamazaki variometer^ for earthquakes at teleseismic 
distances. Most importantly, the experiment proves that georesistivity variations are 
dependable precursors to tectonic earthquakes. 


Field Investigations on Geogesistivity Precursors 

As a sequel to the rock fracturing experiments in laboratories, geophysicists 
attempted to investigate the temporal variations of the resistivities of the crustal rocks 
prior to earthquakes under actual field conditions. They were aware that the success of 
the study depends on the location, depth and dimension of the hypocenters. The length 
of the slipped fault plane is known to vary from about one km in case of magnitude 4 
to about 500 km in case of 8.5 magnitude earthquakes. In spite of uncertainties about 
location of focii, early investigations were carried out at a few sites, over active fault 
zones in Garm regions. Tadzhik republic and San Andreas fault, California, USA. In 
contrast, the Chinese geoscientists monitored earth resistivities at more than- 100 sites 
in the earthquake prone northern China to detect the anomalous focal zones and then 
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to plan intensive measurements to study the precursors of earthquakes. The results of 
these studies have been summarized in the following paragraphs. 


Observations close to the hypocentral regions 

In the Garm district of Tadzhik republic, the apparent resistivities of the crustal 
rocks were measured employing a modified dipole-dipole configuration of electrodes. 
While the electric pulses of about 100 amperes amplitude were passed into the ground 
through a pair of current electrodes^ the potential differences due to this current were 
measured by a pair of electodes at about 6 km from the current dipoles. The ratio of 
the voltage observed by the potential electrodes to the electric current reflects the 
apparent resistivities. The relative or apparent resistivities assuming the value of 30-6- 
67 as reference decreased prior to the occurrence of earthquakes. The maximum 
change of resistivities varying from 7-18 percent was observed to depend on the 
distance between the hypocenter of the earthquake and the center of the dipole unit. It 
was also observed that the magnitude of an impending earthquake is directly 
proportional to the total precursor time. 

j 

Mazzella and Morrison’ measured apparent resistivities on the San Andreas fault, 
California deploying a deep resistivity unit comprising of 85 kW generator, a 
switching rectifier bridge and receivers. The electrodes comprising of aluminium 
plates were buried in pits and the current dipoles were separated by a distance of 1 .5 
km. Out of the potential electrodes, one pair was kept at a distance of about 1 .5 km to 
the north east of current dipoles The 2"“’ pair of potential electrode was kept at a 
distance of about 10 km to the South west of the current dipoles. The 3"^^ pair was kept 
at 1 5 km to the South west, i.e., about 5 km from the 2"‘* pair. The arrays of electrodes 
were arranged along NW-SE direction across the fault plane. A square wave current 
with an amplitude of 200 A and a period of 10s was passed into ground, and the 
potential differences were measured by the receivers. Measurements were carried out 
from January, 1973 to February, 1974. Precursory changes of about 20 percent in the 
resistivity preceded the magnitude 3.8 earthquake on 22"^ June, 1973. These changes 
started since middle of April, 1973. Another activity occurred 8-13 October, but no 
significant changes were observed. 

The above measurements at Park Field have been continued until the present time, 
but no significant variations due to earthquakes have been noticed. This leads to the 
conclusion that any possible changes are below the detection level'^ Considering the 
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resistivity and other precursors, Scholz et.al.'^ proposed dilatancy diffusion model in 
which diffusion of water into the dilatant volume is envisaged and earthquake occurs 
near the maximum stress. According to another model proposed by Myachkin et al 
no diffusion of water is required and in this model, earthquake occurs after the pre- 
anomaly level is reached during a period of falling stress. 

Georesistivity precursors in China 

The Chinese scientists started an elaborate programme of systematic geoelectrical 
investigation establishing a national network of 14 stations in 1966, which increased 
to more than 100 permanent stations in northern China by 1994^ They set up 
additional stations to cover the anomalous areas where impending earthquakes are 
suspected. They carry out repetitive measurements once daily increasing the 
frequency to 6 times per day or continuous if required to study the effects of crustal 
tides. 

In China, georesistivity and geopotential methods have been accepted as most 
effective and powerful techniques for predicting earthquakes. They studied more than 
30 geoelectric precursors to strong earthquakes. The M 7.3 Haicheng e^ent on 4* 
Februaiy, 1975 (the first major earthquake to be successfully predicted), M 6.7 
Xiaojin event on 22"“^ September, 1989 and the M 6.2 Zhang Hei event on lO* 
January, 1998 were successfully predicted to their credit. The M 7.8 Tangshan event 
of 28* July, 1976 was thoroughly investigated, the suspected epicenter was visited 
(Zhao, personal communication), but this disastrous earthquake was not predicted for 
reasons not disclosed yet. However, the wealth of data collected before during and 
after the event has provided insight for refining the interpretation techniques to predict 
other earthquakes. 

In Beijing-Tangshan area within a length of 150 km 8 out of 1 1 stations showed 
persistent decrease, 2 stations showed increase and one station did not show any 
change in the resistivities over a period of 3 years preceding the 1976 Tangshan event. 
The change with greater amplitude appeared earliest around the epicentral zone. The 
anomalies advanced outward with progressively decreasing amplitude from the 
epicenter to the surrounding areas and the maximum amplitude at every station could 
be contoured on a map. According to Zhao et. al.^°, the decrease reflects the process 
of strain accumulation in the area. As shown in Fig. 6, the resistivity anomalies 
appeared late 1973 at Tangshan and Changli (distance-70 km), but after about 3 
months at Baodi (distance-80 km). 
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Fig. 6-^raph of monthly anomalies in resistitivities recorded continuously at six stations around the 
M. 7.8 Tangshan (China) earthquake epicenter. The resistivity decreases were detected 
after middle of 1973 near the epicentral area. The anomalies were observed by the 
end of 1974 at QG and ZXZ stations located 1 10 and 150 Kms from the 
epicenter (after Zhao et. at , 1994). 
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In a new area, the pattern of migration of the anomalies can indicate the location 
of the would-be epicenter Analysing the pattern for many strong earthquakes, Zhao et 
al.* have prepared a travel time curve for the front of resistivity precursors (Fig I Ob) 
that can be used to locate the epicenter using the travel time or duration of onset of 
anomalies observed at three or more stations using the simple technique of 
intersection of circles 

The scientists computed the magnitude of impending earthquake using the 
empirical formulae one each for three parameters namely total duration of the 
anomaly observed at the epicenter, ratio of maximum change in resistivity anomaly to 
pre-anomaly value, and the distance to the farthest station where anomaly was 
detected from epicenter. These formulae have been given elsewhere in this paper. 

The factor for predicting time has been divided into four stages. The long term 
stage depends on the characteristics of the deep crustal conductor mapped by 
magnetotelluric method, which is outside the purview of this paper. The second stage 
or medium term stage starts with appearance of resistivity precursor that develops 
rather slowly. It has been estimated to become about 5 percent over a period of 3 years 
or a rate of 5 x 1 0'^ per day. The speed of change increases with onset of short term 
stage. It was observed more than 3 months before the main Tangshan event when the 
rate increased from 5 x 10'^ to 6.7 x lO”^ per day. The rate increased further in 
imminent stage during the last 15 days and it was 3x10'^ as shown in Fig. 7^°. During 
the imminent stages, the farthest stations exhibited recovery. At the last two stages, 
some stations showed anisotropy, i.e., while the resistivities increased in S-N 
direction, it decreased in E-W direction. 

The critical or imminent stage is associated with some triggering phenomena. One 
of them is the sharp changes or kinks impressed on accelerated S.P. anomaly (Fig. 1). 
These features are accompanied by foreshocks. Under normal conditions, as in the 
case of Bhuj earthquake^', it is difficult to discern the foreshocks. This is rendered 
simpler with awareness created by the geolectrical forerunners. Besides the resistivity 
(Fig. 7z) pulses apparently correlatable with lunar phases show large amplitudes, and 
at this stage, appearance of oil gushes in otherwise old abandoned wells were seen 
near Tangshan. In the final stage, SP may change abruptly before the shock. 

The Chinese scientists believe that monitoring resistivity and self potential are not 
only simple methods for detecting stress-strain within the crust, but also effective 
techniques in the study of earthquake prediction. They caution that negative cases 
have been recorded, which has created trouble with predictions. 
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Georesistivity monitoring in Japan : 


Continuous monitoring of electrical resistivities of rocks has been carried out by 
Yamazaki^ and others in Japan using a very sensitive variometer since 14* May 1968. 
They have passed into ground alternating current having intensity of 100 ma at 250V 
and frequency 67 Hz. The current has been maintained at the odd frequency to avoid 
possible contamination by stray currents from commercial lines in Japan. The 
alternating current has been passed through a pair of outer electrodes in Wenner 
configuration of 4 electrodes while the inner pair of electrodes are used for measuring 
the potential difference. The system has been kept on tuffaceous rocks in a cave at a 
depth of 4.8 m below the ground surface, and electrode separation of 1.6 m has been 
used for measurements. The rate of resistivity change has been recorded continuously 
by a multi-channel plotting type recorder. 



Fig. 8— Record by the Yamazaki variometer for co-seismic and precursory changes associated with a 
M 7.9 16'*' May 1968 earthquake m Japan. The epicenter is located at 800 km from the 
observation point. The step like change nearly coincides with the instant of 
earthquake (Yamazaki, 1977). 
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A rate of resistivity change of the order 10"^ has been observed corresponding to a 
strain value of the order 10'^ and hence with amplification factor of 10^. They have 
observed more than 30 coseismic resistivity indications within 8 years period for 
earthquakes at epicentra! distance vaiy'ing from 60 to 1000 km. Fig. 8 shows a 
coseismic resistivity change for M 7 9 Tokachi-Oki earthquake on 16* May, 1968. 
The rate of change is 0.72 x lO"^ and about 300s are required to complete the change. 
The precursory resistivity change appears to start 4*' preceding the main event. They 
have developed a method to determine the magnitude of the earthquake. They 
observed 12 increases and 18 decreases till 1977 in the trends of resistivities for the 
recorded events. This method can indicate the precursory changes hours in advance at 
teleseismic distance and intensive observation by a number of variometers can help in 
prediction of earthquakes. 

Georesistivity monitoring near Shillong 

The author and his colleagues carried out repetitive resistivity measurements on a 
hill in Shillong, Meghalaya as described under geoelectrica! potentials in this paper. 
The measurements were started on 12-2-1980 with objectives to determine the tidal 
effects on rocks during the total solar eclipse on 16* February, 1980 and to study the 
georesistivity changes as related to earthquakes. The resistivities were measured using 
a potentiometer and accessories manufactured by Geological Survey of India powered 
by battery packs. Besides a Scintrex time domain 2.5 kW I.P unit powered by 
generator was used at the same site repeating observations passing currents at 400, 
600 and 800 volts through a transmitter unit and measuring the potential difference by 
a receiver unit. The separation between the current electrodes along N 40°W direction 
on the slope of the hillock was 380 m with potential electrodes separated by 40m ; and 
the current electrodes and potential electrodes along N20° E (strike direction of rocks) 
were separated by 120m and 20 m respectively. The measurements were continued for 
more than one year although the paper by Nayak et. al?^ describes data upto August, 
1980. 

The resistivities decreased by about 25 percent in some cases and then increased 
to attain pre-anomaly level before earthquakes. Significant earthquakes were observed 
to occur when the resistivity values approached the normal level. As shown in Fig. 4, 
the resistivities decreased gradually and then rapidly from 12,800 Qm to about 10,300 
Qm over a total period of 7 days starting from T' March, 1980. An earthquake M 4.2 
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with its epicenter located approximately 55 km from the site occurred on 9* March, 
1980 at a time when the resistivities were increasing. This trend was repeated for a 
few more events occurring at distance from 40 to 400 km. It was observed that the 
maximaiTi change of the resistivity values is proportional to the hypocenter distance of 
earthquakes (Fig. 10a). 

It was also observed that while the resistivity decreased at Shillong, the values 
i.ncreased for as &TT&y in a bow! like depression at a site about 80 km from shillong for 
a fest earthquake on 29'” .April, 1980. 

Georesistivity monitoring in Maharastra 

Nayak et al}^ took repetitive resistivity measurements at three sites across a 
prominent fault zone separating Gondwana rocks occurring in a depression from 
granite gneisses on a high ground all covered by basalts at Nagpur. The total spread 
lengths of 50m in Wenner configuration of electrodes were laid in the depression and 
on the high ground, whereas the spread length of 320m was laid near the edge of the 
inferred fault. The resistivities observed in July and August, 1991 showed increase of 
54 percent (Fig. 9) in the depression, decrease of 7 percent on the fault edge and no 
change on the high ground. In Killari area of Maharastra, the apparent resistivities 
measured in the epicentral zone of M 6.3 disastrous earthquake of 30* September, 
1993 showed significant changes before the medium magnitude and swarms of 
aftershocks. 

The results of the above investigations suggest that the changes observed by 
repetitive resistivity measurements at sensitive sites can suggest the epicenter and time 
of earthquakes. Since the stations selected on the basis of their sensitivity can detect 
signals due to earthquake upto about 400 km, the number of permanent stations for 
systematic studies can be reduced without sacrificing the advantages of dozens of 
stations as in the case of Chinese techniques. Besides, the false signals can be 
discriminated monitoring the simultaneous increasing and decreasing trends in the 
pairs of configurations as suggested by the author. 

Prediction of Earthquakes 

The prediction of earthquake involves determination of three cardinal parameters 
such as epicenter, magnitude and time of the impending earthquake as accurately as 
possible at the earliest. Although, the procedure has been briefly described under the 
methodology and case histories, it is necessary to synthesize them. 
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Fig. 9 (a>“Temporai variation of magnetic residuals between 1 and 20*^ July, 1991 at a site at 
Nagpur, (b) The apparent resistivities increased by 54% at Gittikhadan on the 
downthrow side of a fault, (d) The resistivities decreased 
by 7% on the flank of the fault zone 

Forecasting epicenter : There are four different geoelectric techniques to forecast 
the epicenters of impending earthquakes (a) Based on the field data of strong 
earthquakes, the Chinese scientists have prepared a travel-time curve (Fig- 10b) that 
shows the relationship between the duration of consistent anomaly at each station and 
its distance from the epicenter^. The epicentral distances from stations showing 
distinct precursors from the travel time curve can help determine the would be 
epicenter with accuracy of 100 km, (b) Contour map prepared on the basis of mean 
decrease as a percentage can also indicate the approximate location of the epicenter 
which will occupy the centre of maximum anomaly, (c) Nayak et al} have shown that 
the distance to the epicenter is related to the ratio of the magnitude of the anomaly to 
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the pl-e-anomaly value. In this method, fewer resistivity stations where the resistivities 
recover before the impending earthquake can indicate the epicenter (Fig. 10a). The 
relative geoelectrical signal densities at a number of stations can be compared to 
determine the epicenter with an accuracy of 100 km'*. 




Fig. 10 (a) The relation between the maximum resistivity change near Shillong vs. distances of 
hypocenters of earthquake (Nayak er al 1983). (b) The travel time curve for the 
resistivity precursor front that can be used to determine the would be epicenter 
in the imminent stage of earthquake. (Zhao et. al 1994). 
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Comparing the techniques, the epicenter can be determined with reasonable 
accuracy combining the Chinese travel time curve and the graph described by Nayak 
et al}. 

Determination of magnitude : The Chinese scientists have developed the 
following formulas on the basis of analysis of data collected for more than 50 
earthquakes to forecast the magnitude of events. 

0.50 + 2.50 logr 
10.80 + 3.41 log A pxipx 

M- -1 .34 + 3 40 log A/h 

Where T is the duration time, A px is the amplitude of the anomaly and p^: 
background resistivity at station near the epicenter, whereas Aw is the maximum 
epicentral distance of the station exhibiting intermediate-term precursor with some 
limitations. The duration of the anomaly may be 2-3 years observed over 200-300 km 
for'earthquake of magnitude M-7. 

Varotsos and Alexopoulos'’ computed magnitude from the relationship between 
the product of relative signal intensity with the distance of the epicenter, but the 
procedure is rather cumbersome. 

Foretelling time of earthquakes 

The long term prediction involving decades are generally based on the calculation 
of return period, analysis of incidence of past earthquake, study of seismic risk zoning 
maps and mapping of the deep crustal conductors. They are useful for creating 
awareness regarding the incidence of earthquake. The intermediate, short and 
imminent stages involving periods of 2-3 years, 2-3 months and 2-3 weeks 
respectively are reflected in the slopes of the graphs depicting the magnitude of 
resistivity values against time. The rates of changes are discernible as in Fig. 7. At the 
short and imminent stages, SP anomalies must be integrated with the resistivity data 
characterised by anisotropy and recovery of anomalies at the farthest stations. At the 
imminent stage, continuous recording of resistivity indications may show pulses of the 
weakened crust correlatable to the frequency of crustal tides (Fig 7z). The amplitudes 
of the pulses may increase distorting the wave shapes and manifesting reversal phases. 
They must be integrated with triggering phenomena such as foreshocks with 
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impressions on the SP anomalies oil gushes in dry wells or water sprouts in tube 
wells. 

The precursory changes observed by Nayak et. al} may correspond to the 
intermediate, short term and recovery phenomena at the stations away from the 
epicentral zone described by the Chinese scientists. The transient changes observed by 
Varotsos and Alexopoulos'* 6-114 h before earthquakes in Greece may be considered 
as a supplementary technique. It may be noted that the Greek scientists announced and 
issued telegraphs predicting earthquakes in Greece and out of 23 only 2 were missed. 
The precursory resistivity changes 0.5-6.8 hr preceding earthquakes by the sensitive 
Yamazaki variometers may also be useful. 

Management of earthquake prediction 

It is known that an earthquake does not kill people, but the buildings which should 
protect people collapse due to faulty constructions and kill people. Hence, the old 
buildings should be strengthened and the new buildings made earthquake resistant in 
all the cities in high seismic risk zones. Since it is not possible to strengthen all the old 
buildings, earthquake prediction programmes must be launched to avert the disasters 
witnessed in 1993 Latur, 1997 Jabalpur, 1999 Chamoli and 2001 Bhuj earthquake. 

It has been shown in the paper that earthquake prediction programmes can be 
initiated employing the inexpensive technique using resistivity potentiometer. The 
seismic zoning map of India shows that Delhi, Kolkata and Mumbai are vulnerable to 
earthquakes. Out of them, case of Mumbai may be considered as an example. 
Earthquakes have been occurring around Koyna to the south of Mumbai and to some 
extent, west of Nasik to north of Mumbai since the nineteen sixties. The disastrous M 
6.3 1993 Latur earthquake occurred to east of Koyna, Hence, the possibility of 
earthquake prone seismic gap in the Western Ghat hills to east of Mumbai cannot be 
ruled out. The presence of hot springs in the area suggest anomalous subsurface 
conditions. Located in Zone-II of the seismic zoning map, Mumbai is also associated 
with prominent high gravity anomalies. The positive isostatic anomaly may indicate, 
deep crustal instability that warrants geoelectrical monitoring on priority basis to 
predict earthquakes. 


Conclusions 


It has been shown in this paper that the gradual accumulation of strain and its 
sudden release to cause earthquakes can be effectively indicated by monitoring the 
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spatial and temporal variation of geoelectric parameters of rocks. Careful analysis of 
the profile of earth resistivities and natural potentials can reveal the epicenter, 
magnitude and time of earthquakes that can not be done by any other method. While 
the Russian and Greek scientists have predicted the instant of earthquake depending 
on geopotential variations, the Chinese scientists have been predicting all the 
parameters integrating georesistivity and geopotential indications at a number of 
stations particularly in the vicinity of selected cities. The Japanese geophysicists also 
claim to have developed method to achieve these objectives employing sensitive 
Yamazaki variometers. 

In India, the author has found that measurements at sensitive sites located on fault 
zones can enhance geoelectric signals due to strain accumulating in the crust within a 
range of 200-400 km from the site and thus a few stations would be needed for 
effective monitoring of parameters. Hence, it is suggested that monitoring of 
geoelectrical parameters near large cities should be undertaken to study the condition 
of crust in order to evaluate the earthquake risk using potentiometers at the earliest 
before it is too late for the next disastrous earthquake. 
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